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SYNOPSIS 

Railway is the most important mode of transportation of 
both men and materials in India. To increase the transportation 
efficiency, heavier, longer and faster trains are needed. The 
track-foundation-soil system is one of the key factors in 

bringing about these changes. The present method of design of 

railway formation is based on empirical basis which gives 
uneconomical design. Moreover, this method does not help 
decision making in respect of suitability of existing tracks 
for heavier and faster trains. 

When trains pass over the track, subgrade gets subjected 
to cyclic stress. If this stress is greater than a certain 
level of stress called threshold stress for soil, cumulative 
permanent deformations are produced leading to ballast 

penetration into the subgrade progressively causing the 
subgrade failure. In some instances, ballast penetration up to 
1 m depth into the subgrade has been observed in the field. 

The present study to predict the behaviour of soils under 
the railway formations subjected to cyclic loading has been 
carried out in three steps viz. (i) mathematical modelling for 
determination of induced stresses in the subgrade based on 



elastic method, 2-D finite element method and 3-D finite 
element method, (ii) geotechnical testing to evolve a suitable 
laboratory testing procedure for formation design and to 
evaluate threshold stress, and (iii) formulation of formation 
design methodology. 

The work contained in the thesis has been organized into 
six chapters. In Chapter One, the topic has been introduced and 
the scope of the present study has been outlined. The available 
literature related to analytical studies for stress 
determination (empirical, finite element and other computer 
oriented approaches) ; geotechnical testing for total stress 
and effective stress studies on both saturated and non 
saturated soils ; track formation design ; and other important 
aspects such as effect of train speed on stresses, 
characterization of track foundation, soil improvement methods, 
field measurements and model testing has been reviewed in 
Chapter Two. 

Mathematical modelling carried out in this study, is 
described in detail in Chapter Three. Three models with 
increasing sophistication viz. ELASTIC. F , 2D8N and 3D20N model 
have been developed for calculating major and minor principal 
stresses in the track structure. 

ELASTIC. F model consists of determining the stresses at 
sleeper-ballast interface by considering the rail as a beam on 
elastic foundation and then using Odemark's method and 
Newmark's solution for calculating stresses in the subgrade. 
2D8N model is a two-dimensional linear, elastic finite element 



model which uses 8-noded rectangular elements. Also, a 
semi-analytical modification to this model is suggested which 
invokes two trials of 2D8N model. 3D20N model is a 
three-dimensional linear, elastic finite element model which 
uses 20-noded brick element and models rail as one-dimensional 
beam element. Performance of this model is very similar to 
'MULTA' model [Selig et al. (1979)] which is like a closed form 
solution. 

Parametric studies using all the three models have been 
carried out to evaluate the relative importance of each track 
component. It is observed that depth of formation and type of 
subgrade are the most important factors. 

Results predicted by all the three models have been 
compared with predicted results by other models ; measured 
results in the field by ORE (Office of Research and 
Experiments, International Union of Railways) and Indian 
Railways ; and full scale model test results by ORE . It has 
been observed that all the three models predict quite 
satisfactory results. Comparisons have been also made with 
values predicted by semi-empirical methods and it is shown that 
these methods overpredict by as much as 2 to 4 times the 
measured results. 

It is shown that vertical and horizontal stresses 
predicted by elastic analysis compare well with those obtained 
by elasto-plastic analysis which takes into consideration 'no 
tension' condition in granular media. A non-dimensional 
parameter 'm' has been introduced as a measure of ballast. 



subballast and subgrade characteristics (i.e the formation 
quality) . 

Chapter four deals with the geotechnical testing carried 
out on alluvial silt for evaluation of strength, 
stress-strain-pore pressure behaviour during static and cyclic 
loading. Total six series of tests viz. total stress static 
shear tests, unconfined cyclic tests, unconfined cyclic 
softened tests, effective stress static shear tests, cyclic 
undrained tests on quasi saturated samples with pore pressure 
measurement at high confining pressure and cyclic undrained 
tests on back saturated samples with pore pressure measurement 
at low cell pressure simulating the track loading conditions, 
have been carried out. 

Relationship between incremental plastic strain, Se^ after 
100 cycles on unconfined compacted samples and corresponding 
increase in incremental water content on softening has been 
indicated. The explanation of mudpumping problems encountered 
by Indian railways under the tracks in alluvial soil deposits 
has been well brought out by observing the significant 
degradations in both stiffness and strength following the 
transient loading and softening. 

A simple procedure to ensure quasi saturated response for 
partially saturated sample is given. It has been shown that 
stress - strain, pore pressure - strain response and effective 
stress paths for compacted samples follow non-normal ized 
behaviour with respect to consolidation stress. It is observed 
that evaluation of failure envelope by Marsal's procedure gives 



satisfactory results. 

Stress - strain, pore pressure - strain , stress path and 
other results for IITK Campus silt show remarkable similarities 
with group IV soils [Cruz et al. (1985) ] and with Canyon dam 
clay [Casagrande and Hirschfeld (1960) ; Casagrande and Poulos 

(1964)] in all aspects of behaviour of compacted samples. 
Effective stress data available for Campus silt and Canyon dam 
clay on compacted samples has been examined in the general 
framework of critical state model. It has been brought out that 
compacted soil samples behave like overconsolidated soils with 
different values of past maximum pressure (produced due to 
compaction, depending upon moulding water content and dry 
density for a given compaction energy) . It has been observed 
that compacted soils, when tested at a degree of saturation 
above a certain value which ensures quasi-saturation, behave 
like over consolidated soils up to effective confining stresses 
below that needed for complete saturation and like normally 
consolidated when tested at effective confining stresses 
greater than that needed for full saturation. 

Guidelines given by Yudhbir and Wood (1989) have been 
followed in drawing Virgin consolidation line, Critical state 
line, Recompression lines, Tension cut off etc. It has been 
observed that these lines show reasonable fit with actual 
experimental data. 

The generalized relationship, initially proposed by 

Leroueil et al. (1992) between (cj-w , )/I and s has been 

opt p u 

modified. This relationship may be used for estimating 
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undrained strength of compacted clays at OMC or wet of OMC. 

A generalized relationship to relate cumulative plastic 
strain during cyclic loading with number of cycles has been 
proposed and the values of coefficients reguired in this 
relationship for different types of soils have been evaluated. 
It is observed that end points of cyclic load tests on 
compacted samples would lie on Hvorslev surface as observed by 
Hyde and Ward (1984) for normally and overconsolidated soils. 

Procedure to evaluate threshold stress on the basis of 
cumulative plastic strain data is suggested. Rubin et al. 
(1970) also suggested a similar procedure on the basis of 
drained cyclic tests. Threshold stress has also been evaluated 
from the equilibrium line concept suggested by Sangrey (1968) . 

A generalized relationship between threshold stress ratio, 
R,_ 0 and plasticity index, I has been developed for a variety 
of soils (compacted, stiff, normally consolidated etc.) at low 
and relatively high effective confining stresses. Relationship 
of the form (o^-tTg)^ = a + b a ' has also been given for a 
variety of soil types at low effective stress ranges. 

In Chapter Five, a rational design methodology based on 
stresses induced and threshold stress evaluated from above 
study has been proposed. A detailed step by step procedure 
(flow chart) has been outlined and a simplified method for 
routine design of formation depth has been recommended. 

The conclusions arrived at on the basis of present study 
and the recommendations for the further studies are given in 
Chapter Six. 



CHAPTER 1 


INTRODUCTION 


1.1 GENERAL 

Railway is the most important mode of transport of both 
men and materials in India. With the economic development and 
growth of industries, need to increase the transportation 
throughput is being increasingly felt. Transportation 
throughput can be increased by running longer, heavier and 
faster trains besides other measures. The capacity of 
track-foundation— soil system is one of the key factors in 
bringing about these changes. 

At present, formation thickness under the railway track is 
decided on empirical basis. This method on one hand, 
overpredicts the thickness where subsoil condition is good and 
is, thus, uneconomical while on the other hand, underpredicts 
the formation thickness where subsoil condition is poor 
requiring frequent and expensive maintenance measures and 
obstruction to the flow of traffic. Moreover, this method does 
not help decision making for heavy axle loads and high speeds. 
Thus, in order to increase the speed and the load carrying 
capacity of trains, it is imperative that a new design 
criterion based on a rational approach be formulated. This 
study would also help to establish the suitability of the 



existing formations with respect to their capacity to carry 
heavy and high speed traffic. 

1 . 2 PRESENT PROBLEM 

Railway formations in India were constructed till very 
recently on procedures based on the experience of the engineers 
involved. A slope of 2:1 for banks and steeper slopes for 
cuttings were considered to be safe. The ballast depth and 
sleepers provided for the track foundation were generally 
sufficient because of the lighter axle loads and less frequency 
and speed of trains. 

However, as the traffic developed, heavier axle loads, 
longer and faster trains and increased frequency created 

problems of track maintenance. This necessitated the 

A 

development of better track structure viz. rails, sleepers, 
sleeper density and ballast cushion. However, despite these 
attempts several large stretches of the track could be 
distinguished where the subgrade soil was responsible for 
instability of track structure / foundation. Moreover, 
development of better track structure was not based on any 
rational design. 

On Indian railways, the need for treating such locations 
having subgrade failure was felt when high speed Rajdhani 
express trains were planned to be introduced. Simultaneously, 
on heavily worked routes of Eastern railways, selected remedial 
measures were tried. Many of these remedial measures were 
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picked up from procedures and practices evolved for the new 
Tokaido line in Japan. In the following years, behaviour of 
treated locations indicated the intricacy of the subject and 
many of the formation treatment methods which seemed to be very 
promising in the beginning, did not solve the problem (Table 
1 . 1 ) . 

In general, causes of the subgrade failure mentioned above 
can be attributed to the following causes : 

a. Poor construction methods (especially lack of proper 
compaction control) . 

b. Lack of proper drainage of track formation. 

c. Empirical nature of design methodology without any 
systematic regard for the effect of cyclic loading on the 
behaviour of track-foundation-soil system. 

Figure 1.1 illustrates the broad outline of subgrade 
failure mechanism. When trains pass over the track, subgrade is 
subjected to a certain cyclic stress. If this stress is greater 
than a particular stress level called threshold stress for the 
soil, ballast penetrates into the foundation soil progressively 
causing the subgrade failure. Penetrated ballast tends to store 
water percolating through the track surface. Under the train 
loads, mud pumping phenomenon is observed. Mud pumping of track 
is a phenomenon caused when roadbed soil is softened by the 
rainfall or groundwater; the mud climbs up through the void 
inside the ballast and pumped out on the ballast surface by a 
passing train load. This displaces the clay and fines present 
in subgrade soil and completely disrupts the drainage system 
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Table 1 . 1 Status of formation Treatment Methods 


No. 

Methods of formation Treatment 

Status 

1 

Lime Piles 

* 

IE 

2 

Vinyl Drains 

IE 

3 

Sand Drains 

IE 

4 

Blanket Moorum 

PE 

5 

Polythene sheet 

IE 

6 

Cationic Bitumen Emulsion 

IE 

7 

Cement Pressure Grouting 

PE 

(Limited Life) 

8 

Maxphalt - Crete 

IE 

9 

Ballast Piles 

UT 

10 

Partial Blanket 

UT 

11 

Lime slurry pressure Injection 

PE 

12 

Geotextile 

IE 

13 

Pouring Sand on Ballast 
(Protzeller 7 s Method) 

PE 

14 

Blanket - Coarse grained 

FE 

15 

Later ite Blocks 

UT 

' 

Method of Execution 


1 

Blanketing under traffic 

PE 

2 

Relieving Girder 

UT 


* IE = Ineffective 

FE = Fully Effective 
PE = Partly Effective 
UT = Under Trial 










Track structu 




causing further track deterioration. This results in increase 
in costly and difficult maintenance cycles, making track 
operation uneconomical. If proper and timely maintenance is not 
carried out, it may also lead to serious accidents. 

The present study aims to evolve a rational criterion for 
design and construction of railway formations by studying the 
effect of cyclic loading on subgrade soil. 


1.3 SCOPE OF THE PRESENT INVESTIGATION 

The proposed study has been carried out in the following 
steps . 

1. Mathematical modelling of railway track to determine the 
stresses 

These studies are carried out in order to estimate the 
pressure on top of formation due to the faster and heavier 
train loads using Finite Element method. Parametric studies are 
carried out to evaluate the relative importance of each 
component of the track. Design charts for determination of 
induced stresses in the subgrade have been prepared. 

2. Geotechnical testing for evaluation of strength, stress - 
strain - pore pressure behaviour 

This is done in order to evolve a suitable laboratory 
testing program to simulate track loading conditions in the 
subgrade soil and to develop procedures for estimation of 



7 


threshold stress. General guidelines for evaluation of 
threshold stress for a variety of soils are suggested. 

3. Formulation of design methodology 

Based on the analytical and experimental studies, a design 
methodology for track formation has been recommended. 

The details of the present investigation have been 
presented in Six chapters as outlined below. 

Chapter I presents the introduction of present problem 
indicating the need of present study and broad outline of 
investigations carried out. 

Available literature on the analytical studies for stress 
determination, geotechnical testing related to track loading 
conditions, track formation design and other important aspects 
related to present problem has been reviewed in Chapter II. 

Chapter III deals with mathematical modelling of the track 
structure and the results obtained have been analysed. 

Details related to geotechnical testing program, test 
results and interpretation under both static and cyclic loading 
conditions have been discussed in Chapter IV. 

In chapter V , formulation of design methodology has been 
presented. 

The conclusions arrived at and the recommendations for the 
future studies are given in Chapter VI. 



CHAPTER 2 


LITERATURE REVIEW 


2.1 INTRODUCTION 

In order to predict the behaviour of soils under the 
railway formations, it is imperative that the stresses inside 
the track structure should be estimated by proper mathematical 
modelling and the loading condition should be simulated in 
laboratory to obtain the threshold stress. 

The trends observed in literature related to mathematical 
modelling can be broadly divided into models which incorporate 
the exact solution for multilayer elastic system given by 
Burmister (1945) and models based on finite element analysis 
(elastic and elasto-plastic) . However, all models give the 
stresses assuming equivalent static load on top of rail. This 
equivalent static load is obtained by augmenting the wheel load 
by suitable speed factor. 

The simulation of track loading condition in laboratory to 
determine the threshold stress for soil has been done by Rubin 
et al. (1970) for London clay and following this approach, 
Janin et al. (1983a) have also determined the threshold stress 
for Keuper marl and Boulder clay. 

The literature review in this chapter has been carried out 
under the following heads : 



a. Literature related to analytical studies for stress 
determination . 

b. Literature related to geotechnical testing for evaluation of 
strength, stress - strain - pore pressure behaviour. 

c. Literature related to track formation design. 

d. Literature related to other important aspects such as speed 
effect, characterization of elastic foundation, soil 
improvement methods, and field measurements and model testing. 

2.2 LITERATURE RELATED TO ANALYTICAL STUDIES FOR STRESS 
DETERMINATION 

Stresses inside the track structure were earlier 
determined by empirical and semi empirical approaches. Finite 
element method and other computer oriented approaches based on 
closed form solution of multilayer system given by Burmister 
(1945) were increasingly used for stress determination after 
1978. This section has been classified under empirical 
approaches and finite element and other computer oriented 
approaches. 

2.2.1 EMPIRICAL APPROACHES 

1. Zimmermann's Method for Stresses in Rails 

Based on assumption of rail supported on infinite number 
of elastic bearings, Zimmermann gave expressions for 
deflection, moment and stresses on rail as given below : 
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y r = Qtj(€)/2 

(2.1) 

M = Q L b U(^)/4 

(2.2) 

ff rl - fi / w l 

(2.3) 

a r2 - M/W 2 

(2.4) 

S - F / (2 S) , 

(2.5) 

• 

b C 

(2.6) 


v = deflection of rail , 
r 

Q = wheel load , 

M = total moment on rail , 

x = co-ordinate of point where deflection is needed, 
C = formation coefficient, kg/cm , 

€ = x/Lfc , 

b = bearing factor , 

S = spacing between sleepers , 

F = bearing area per sleeper , 

= base length of track , 

E = Young's modulus of elasticity of rail , 

I = moment of inertia per rail , 

= moment of resistance of rail foot , 

W 2 = moment of resistance of rail head , 



XX 


a . = stresses on rail foot , and 
rl ' 

cr r2 = stresses on rail head. 

Charts to pick coefficients v(^) and h(£) are given in RDSO(B) . 


2 . Method of Odemark for Vertical Stress Calculation in 
Multilayer Systems 

Effect of layering is to reduce the stresses at a given 
depth. In this method, effect of layering is taken into 
consideration by increasing the depth based on elastic modulus 
of different layers and then, calculating the stresses by 
Boussinesq's method at this new depth [RDSO (B) ] . 

Referring to Fig. 2.1 , the expressions for modified depths are 
given below : 


(2.7) 

( 2 . 8 ) 

where 

*k "k , , 

dfo, = modified depths of ballast and subballast 

respectively, 

dfaf d sb = original depths of ballast and subballast 
respectively, and 

Efa, Esk' E s = Young's modulii of elasticity of ballast, 

subballast and subgrade respectively. 


d b = 0 • 9 d b 3 / V E s 


d sb " °* 9 d sb 3 / E sb/ £ s 




Layered system Homogeneous system 


Fig. 2.1 METHOD DF DDEMARK 
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3. Burmister's Multilayer Approach 

Selig et al. (1979) have developed an approach that uses 
Burmister's multilayer theory [Burmister (1945)] for 
calculating stresses on ballast and subgrade and a structural 
model for rail-tie interaction. Contact between tie and ballast 
is represented by a series of circular areas with uniform 
pressure. 

4. Japanese Approach 

Based on elastically supported beam model, pressure, 
bending moment and deflection in the rail are worked out [RDSO 
(1987) ] . 

Empirical stress distribution over the ballast layer has 
been proposed from which average subgrade pressure can be 
worked out. Maximum subgrade pressure can be found out as : 


vt 


vmean 


58 


b 1 


10 + d 


1.35 


b + 2(d b - 15) 


) ( 


1 + 2(d b - 15) 


(2.9) 


where 


b,l = 


d b = 


vt 


O' 


vmean 


width and length of sleeper respectively, 
ballast depth, 

maximum subgrade pressure, and 
average subgrade pressure. 
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5. Formation Pressure based on Technical Paper No. 245 

On Indian railways, the maximum intensity of pressure on 
subgrade is calculated as given below [Kapoor et al. (1967a)]. 


P 

s 


= Q 


64 E I 
r r 


S 


( 2 . 10 ) 


cr = 
Vt 


71 d, L 
b 


( 2 . 11 ) 


where 


P g = load on individual sleeper, 

L = effective bearing length of sleeper, and 
k = modulus of rail road reaction. 


6. Formation Pressure Intensity by C.W. Clarke (Technical 
monograph number 12) 

Both average and maximum pressure intensity can be 
obtained by the following relationship given by C. W. Clarke 
[Kapoor et al. (1967a)]. 


vmean 


0.39 ZLV S 

X, 0.58 b L d, 
1 b 


10 


vt d. Area of sleeper 


( 2 . 12 ) 


(2.13) 
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where ZLV = Zimmermann Wheel load, 

= distance from load on rail to point of zero 
bending moment on rail, and 
P g = to be found from Eq. (2.10). 

7. Talbot's method for American Standard Gauge Track 

A. N. Talbot has proposed following expression for 
pressure intensity [Kapoor et al. (1967a)]. 

°vt 16.8 

= (2.14) 

cr .1.25 

S d b 

where <r = maximum pressure at the bottom of sleeper. 

8. Pressure Intensities Calculated by German Railways 

Kapoor et al . (1967a) have reported the following 

expression for maximum subgrade pressure given by German 
railways : 


1.5 P 


cr = 
V 


3(l-G)+b d, tan(a ) 

c J JJ ± 


(2.15) 


where G = track gauge, and 

a^j^ = angle of distribution of pressure in formation. 



However, it has been observed that the modern trend for 
determination of formation stresses is based on the use of 
finite element and other computer oriented approaches. 

2.2.2 Finite element and other computer oriented approaches 

The various mathematical models that have been used for 
investigation of track design parameters are given below. 

1. So model 

It consists of following four different routines [So 
(1978) ] . 

a. Finite element vertical track model 

It is a two dimensional (2-D) finite element model 
tfhich uses 1-D members to simulate rails and tie foundation 
stiffness. It predicts vertical rail-tie reactions from given 
/ertical wheel - rail loading. 

b. Finite element lateral track model 

It is similar to the above model but it simulates 
:rack structure behaviour under lateral loading and predicts 
Lateral tie loads. 

c. Finite element tie model 

It is a 2-D model with similar formulation and 
:apabilities as above models but it generates vertical and 



lateral tie-ballast loading from vertical and lateral rail-tie 
loading output from the above models. 

d. Multilayer elastic model 

Tie-ballast loading from above is used as input. This 
model is based on Burmister's solution [Burmister (1945)] for 
multilayer elastic system. 

When all the above four models are used in combination, 
they give stresses at any point in the track structure, 
provided compatibility of displacements and forces is satisfied 
at boundary interfaces. 

2. MULTA model 

The multi-layer track analysis (MULTA) model is a 
combination of two routines - Burmister and LAC [Selig et al. 
(1979) ] . 

Loads and Combination (LAC) routine is a matrix structural 
analysis model which solves for the tie-ballast reactions using 
the method of consistent deformations. Rails are assumed to be 
supported by certain number of ties and each with certain 
number of divisions. The tie-ballast contact pressures are 
assumed to be transmitted through these divisions. Reactions 
between rails and ties and displacements at intersection of 
rails and ties are also assumed as unknowns. 

Burmister routine is based on Burmister's multilayer 
elastic theory [Burmister (1945)] and above tie-ballast 



reactions are taken as input and stresses at required points 
are found out. 

3 . PSA model 

It consists of two separate algorithms viz. PSA (Prismatic 
Solid Analysis) code and Wheels code [Selig et al. (1979)]. 

Wheels code is a matrix structural algorithm and is 
similar to LAC routine. 

PSA code generates the foundation stresses and 
displacement coefficients. It performs the analysis of a class 
of 3-dimensional (3-D) problems whose solutions can be 
approximated by analysis of periodically loaded prismatic 
solids. The 3-D solution is approximated as Fourier series 
along the track direction. The coefficients in the series are 
obtained from 2-D finite element analysis. Summation of this 
series yields the final solution. The output consists of 
displacements, strains and stresses at any desired point in the 
track structure. 

4. ILLITRACK model 

It simulates the track system with two pseudo plane strain 
2-D finite element analysis [Selig et al. (1979)]. The 
longitudinal 2-D analysis is followed by transverse 2-D 
analysis using the results from the previous longitudinal 
analysis as input. It employs following two types of finite 


elements . 



a. Rectangular plane strain elements for ballast. 


subballast and subgrade. 

b. Beam-Spring elements to represent rail tie sub-system. 


5. GEOTRACK model 

GEOTRACK [Chang et al. (1980)] is a 3-D, multilayer model 
for the elastic response of railroad track and it considers 
stress-dependent material properties and separation of tie and 
ballast. The computer program 'GEOTRACK' was developed starting 
with the program MULTA and deficiencies noted in 'MULTA' were 
removed. The changes pertain to the following aspects. 

a. Excessively high stresses previously calculated beneath 
the center of each tie were corrected. 

b. Vertical tensile stresses between ties and ballast were 
avoided . 

c. Effects of two axles could be superimposed. 

6. Finite element studies by Desai and Siriwardane (1982) 

Desai and Siriwardane (1982) have presented 3 finite 
element formulations based on one-, two-, and three-dimensional 
idealization for analysis of rail road track support 
structures. 3-D formulation incorporates nonlinear constitutive 
behaviour by using four different elasto-plastic laws viz. 
Variable modulii, Drucker-Prager model. Cap model and Critical 
state model. It uses thin interface elements for simulation of 
behaviour at junctions. 
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a. 1-D model : It assumes rail supported by equivalent 
springs as a beam on nonlinear foundation. Rail is modelled as 
1-D beam with six degrees of freedom system. 

b. 2-D model : It involves plane strain idealization and 
uses 8-noded isoparametric elements for ties, ballast, 
subballast and soil. 

c. 3-D model : Hexahedral elements were used for various 
track components. Bottom boundary was assumed to be rough and 
lateral boundaries as smooth. Boundary distance was fixed by 
parametric studies. A comparison of results was presented using 
PSA, MULTA, ILLITRACK and three dimensional finite element 
model with different boundary distances. However, it is 
observed that results obtained by 3-D finite element are lesser 
than results obtained from MULTA model. 

CPU time comparison has also been given which shows CPU 
time of 1 min. and 44 sec. for linear 3-D analysis and 15 min. 
and 36 sec. for nonlinear analysis on IBM 3032 computer system 
for the mesh adopted therein. 

7. Studies by ORE (Office for Research and Experiments, 
International Union of Railways, Utrecht) 

Janin et al. (1983b) have carried out track analysis by 
two computer programs as given below. 



a. BISTRO model 


: It is an elastic model and is based on 


Burmister's multilayer theory. 

b. ROSALIE model : Both elastic and elasto-plastic , 3-D 
finite element models were developed. They employ hexahedral 
elements with 20 nodes. 

Track analysis using both programs was carried out for 
single sleeper box at Vienna and Derby and for the complete 
panel of track at Derby. Results were compared with field tests 
and single sleeper box tests. 


2.3 LITERATURE RELATED TO GEOTECHNICAL TESTING FOR EVALUATION 
OF STRENGTH, STRESS-STRAIN-PORE PRESSURE BEHAVIOUR 

In order to predict the behaviour of soils under the 
railway formations and to evaluate a rational design procedure 
for railway formations, it is important to study the effect of 
cyclic loading on samples in laboratory simulating the track 
loading conditions. 

This section has been classified into total stress 
/fatigue studies and effective stress studies. 

2.3.1 TOTAL STRESS/FATIGUE STUDIES 

1. Bishop and Henkel (1953) observed the effects of 
applying the stress pulses to undrained samples of Weald clay. 
They observed that residual pore pressures were left after the 
removal of load and that the magnitude and sign of these pore 



pressures were dependent on the consolidation history of the 
clay and on the magnitude of the transient stress. 

For a normally consolidated clay which develops +ve pore 
pressures in shear, they have observed that upon removal of the 
shear stress, a +ve residual pore pressure persists. The 
reapplication of the shearing stress causes the soil to develop 
additional +ve pore pressures, apparently increasing the pore 
pressure to a greater magnitude than would have developed had 
the first cycle of load been carried to failure. 

When a heavily overconsolidated sample was subjected to a 
shear stress increasingly negative pore pressures were 
developed after the development of a small positive pore 
pressure at very low levels of stress and strain. A substantial 
negative pore pressure remained after unloading. 

2. Sangrey (1968) reviewed the existing literature and 
noted that up to 100,000 cycles, the specimen deformation 
depended only on total number of cycles and was independent of 
frequency of cycles (tested for 1 to 20 cycles/min. frequency 
range). For thixotropic clays, frequency affected the 
deformation. 

3. Janin et al. (1980) performed cyclic triaxial tests of 
3 soils namely silt, silty sand and uniform fine sand. They 
made the following observations : 


a. For tests conducted on 10 cm diameter and 20 cm high 



samples in frequency range of 1 to 4 0 hz, no resonance 
phenomenon was observed. Hence they carried out all other tests 
with 5 Hz frequency. 

b. In fatigue tests, the failure criterion adopted was 
plastic deformation exceeding the failure deformation values 
determined in the static triaxial tests. The K-value was 
defined as pulsating load causing failure to the static failure 
strength. All tests were carried out with drained conditions. 

c. It was found that K-value varied with soil type as 
given below. 

K = 0.60 for silt, 

= 0.28 for silty sand, and 

= 0.30 for uniform fine sand, 
for cr 3 = 15 kPa. 

d. Results of fatigue tests at different lateral stresses 
showed a distinct increase of plastic deformation with 
decreasing lateral stress. 

e. More the initial density, lesser the plastic 
deformation in fatigue tests. 

f. Water content affected K-value and deformation to a 
considerable extent in silt as compared to silty sand and 
uniform sand. 

4. Janin et al. (1983a) carried out laboratory tests to 
determine the properties of sandy silt and silty clay used in 



full scale model laboratory tests at Vienna arsenal (Sec. 
2.5.3.) and from these tests, they made the following 
observations . 

a. A difference in water content has more pronounced 
effect on engineering properties of silt than of the clay due 
to low plasticity. 

b. Smooth stress strain curves were obtained for silty 
clay. But sandy silt showed brittle behaviour at low water 
content and showed behaviour of strain hardening materials at 
high water content. 

c. Sandy silt showed extremely high friction angle (45 
degrees) for undisturbed samples due to pronounced dilatancy at 
low confining stresses (<x 3 = 35 kPa) . However, a lower value of 
34 degrees was obtained for compacted samples. 

d. The stress-strain relationships of the samples loaded 
monotonically to failure after cyclic loading (5 x 10 load 
cycles at different cyclic stress levels) were similar for both 
soils. A more brittle behaviour and a lower failure strain were 
observed. 

e. Secant modulus decreases with increasing confining 
stress level. 

2.3.2 EFFECTIVE STRESS STUDIES 

The literature related to effective stress studies has 
been reviewed under the following sections : 



a. Saturated soils 


b. Non saturated soils 

2. 3. 2.1 Saturated soils 

1. Sangrey (1968) conducted the series of experiments on 
saturated clays and made the following observations : 

a. For normally consolidated samples, there exists a 
critical level of repeated stress (CLRS) such that when cycled 
stress level is below CLRS, plastic deformation is terminating 
(i.e. equilibrium condition is reached) and when cycled stress 
level is above CLRS, plastic deformation is non-terminating. 

b. Figure 2.2 shows the typical comparison of 
stress-strain curves of single loading and repeated loading for 
normally consolidated Newfield clay. Repeated loading curve has 
been defined by tests conducted at different cyclic stress 
levels. 

c. Figure 2.3 shows peak strength relation with level of 
cyclic stress. 

d. Tests were conducted on normally consolidated samples, 
anisotropically consolidated samples and overconsolidated 
samples with overconsolidation ratio, OCR =2, 4, 6. Behaviour 
for these samples remains similar to normally consolidated clay 
in general except that as OCR increases, clay exhibits lesser 
residual pore pressure development at a given stress level. 
Typical curve for OCR = 6 is given in Fig. 2.4. 

e. Procedure to evaluate critical level of repeated stress 


“PME^ATERTPRESSURE “ a'u PS I STRESS' Cor-cs) PS I 



1 90. Consol.; 3 14.2 psi, OCR - 6 } 
Single LoadIvs. Rep eaTed~ Load - "| 


single ‘‘{load 


/r^peat^d lokd 


load 


I 

Lie— 


jL jSTRA^N 1 

I 


Figure 2.4 j Compiriso^ of jstreEi-Str&in Curves of 

I i Single Li>adink and! Repeated Loading 




(or threshold stress) based on movement of effective stress 
path due to development of pore pressure during cyclic loading 
was suggested. 

2. Rubin et al. (1970) investigated clay behaviour under 
dynamic and repeated triaxial loading. Standard repeated 
triaxial load test is defined as a test at a constant ambient 
pressure of 35 KPa with cyclic frequency of 30 cycles per 
minute. Failure criterion was assumed as 10 percent cumulative 
deformation produced in less than 1000 load cycles. Threshold 
stress is defined as a maximum stress that can be applied to 
the sample that does not cause cumulative strain greater than 
10 percent in 1000 cycles. 

Based on studies related to development of plastic strain 
with number of cycles at different cyclic stress ratio, a 
relationship between cumulative plastic strain and number of 
cycles has been proposed. 

It has been observed that as cyclic stress increases, the 
nature of relationship between plastic strain development and 
number of load cycles changes from terminating to the 
non-terminating one. Based on this observation and the above 
mentioned failure criterion, Rubin et al. (1970) have 
calculated the threshold stress in drained tests. 

No relationship is apparent between the water content 
changes and the applied stress. 

Threshold stress value is not a basic material parameter 



but is dependent on mean effective stress, wave shape, 
frequency and previous loading history. 

3. Gupta (1973) carried out static and repeated triaxial 
tests on saturated normally consolidated Campus silt. Confining 
pressure was found to have dominant effect on nature of 
stress-strain curve. For repeated load tests, it is found that 
for 25 percent and 55 percent stress level, pore pressure 
increases during loading and decreases with unloading. For 75 
percent stress level, the response is reverse. It is observed 
that for the silty soil on which tests were carried out, 
equilibrium line exists. 

4. Wilson and Greenwood (1974) carried out cyclic tests 
with pore pressure measurements on lacustrine clay at Hamilton, 
Ontario. They presented the data for plastic strain development 
and pore pressure increment with number of load cycles. 

5. Schjetne (1976) has carried out stress controlled 
triaxial tests and simple shear tests on clay with different 
OCR and concluded that c' and <p' remain unaffected by cyclic 
loading where c' and <f>' are the effective cohesion and angle of 
friction of clay. The undrained shear strength after cyclic 
loading is a function of cyclic shear strain reached and the 
number of cycles. 
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Based on undrained triaxial tests on remoulded clay, 
Hyde and Brown (1976) proposed the following relationship . 

log e = C - D log T (2.16) 

where t = strain rate per unit time, 

C = log of strain rate at unit time, 

D = strain rate decay constant, and 
T = time in secs. 

D was found to be constant for samples with the same 
stress history. 

7 . Based on results from laboratory testing of Drummen 
clay under cyclic loading, Andersen et al. (1976) have 
concluded that for normally consolidated clays, drainage will 
reduce cyclic displacements and improve stability while for 
overconsolidated clays, drainage will increase cyclic 
displacements and will reduce the stability. 

8. Based on undrained static and cyclic triaxial tests on 
normally consolidated plastic clays, Yudhbir and Rehman (1977) 
made the following observations. 

a. Clay with lower plasticity index (I p ) has a lower 
value of critical stress level than clays with higher I p . 

b. x cyc / s u (applied cyclic stress, x cyc as a fraction 
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of undrained shear strength, s^) increases with the increase in 
confining stresses. 

c. Pore pressure set up at equilibrium state increases 

with cycled stress level and is higher for low I clay than for 

P 

high 1^ clay. 

9. Sangrey et al. (1978) have presented the relationship 
between threshold stress ratio and k/v^ where /c is the slope of 
swelling line in triaxial consolidation test and v^ is the 
specific volume of sample. The term jc/v ^ is the measure of soil 
compressibility. From this relation, threshold stress ratio for 
variety of soils can be evaluated. However, it must be noted 
that evaluation of k poses certain problems. 

10. Andersen et al. (1980) carried out cyclic triaxial and 
simple shear tests on Drummen clay for OCR from 1 to 50 and 
concluded that the behaviour of clay during undrained cyclic 
loading could be explained in terms of effective stresses. Some 
of the observations made by them are as follows. 

a. A normally consolidated clay is more resistant to 
undrained cyclic loading than an overconsolidated clay. 

b. Results from strain controlled cyclic loading do not 
coincide with the stress controlled cyclic loading except when 
interpreted in terms of accumulated cyclic shear strain. 
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that evaluation of k poses certain problems. 

10. Andersen et al . (1980) carried out cyclic triaxial and 

simple shear tests on Drummen clay for OCR from 1 to 50 and 
concluded that the behaviour of clay during undrained cyclic 
loading could be explained in terms of effective stresses. Some 
of the observations made by them are as follows. 

a. A normally consolidated clay is more resistant to 
undrained cyclic loading than an overconsolidated clay. 

b. Results from strain controlled cyclic loading do not 
coincide with the stress controlled cyclic loading except when 
interpreted in terms of accumulated cyclic shear strain. 



11. Wood (1982) has carried out review on laboratory 
investigations of the behaviour of soils under cyclic loading 
3nd following are some of his findings. 

Effective stresses and pore pressures provide the key 
controlling the behaviour of soils under cyclic loading. 

Test procedures need to be such that pore pressures can be 
reliably monitored. This leads to prefer slow cyclic tests to 
fast ones on the grounds that pore pressure measurements in the 
former are more likely to reflect the actual intentions of the 
soil . 

Only first few cycles of cyclic loading are important to 
make observations on both changes in pore pressure and the 
stiffness . 


12. Hyde and Ward (1985) have shown that permanent 
pressures built up during cyclic loading of normally 
consolidated Keuper marl can be expressed by the following 
relation. 

log (u/P e ') = log a - 0 log N (2.17) 


where 


a = pore pressure rate at n = 1, 

/3 = gradient of log (u/P ) vs. log N plot, 

u = rate of change of pore pressure with N, and 
P ' = Hvorslev's equivalent consolidation pressure. 


They also describe that though this relation is similar to 



Eq. (2.16) but there is no simple relationship between strain 
and pore pressure development. 

They have observed that Hvorslev surface forms the upper 
boundary for cyclic effective stress tests for failure. 

13. Shew et al. (1987) have conducted a series of cyclic 
triaxial tests on Kaoline clay under 3 triaxial stress paths 
viz. cyclic compression, cyclic extension and cyclic 
compression extension and concluded that critical stress level 
for repeated loading for cyclic compression test is larger than 
the cyclic extension tests and is much larger than cyclic 
extension compression tests. 

14. Based on undrained triaxial and direct simple shear 
tests on Drummen clay, Andersen et al. (1988) concluded that 
cyclic properties of clay depended upon both average shear 
stress and type of test (e.g. triaxial or direct simple shear) . 

2. 3. 2. 2 NON SATURATED SOILS 

1. Bishop et al. (I960) gave the following equation for 
shear strength evaluation of partially saturated clay. 

x = c'+(cr - u a )tan <p' + X (u a ~u w ) tan<£' (2.18) 

where x = Undrained shear strength of partially saturated 


soil. 



X 


= a parameter that varies with the degree of 
saturation of sample, 
u = the pore air pressure, and 

CL 

u^ = the pore water pressure. 

2. Lambe (1960) gave mechanistic picture of shear strength 
of clay and proposed alternative relationship for shear 
strength prediction of partially saturated clay based on 
internal attraction and repulsion of clay particles. However, 
this relationship remained mostly academic in nature due to 
difficulty in measuring the internal attraction and repulsion. 

3. Casagrande and Hirschfeld (1960) and Casagrande and 
Poulos (1964) have done extensive studies on stress 
deformation, strength and pore pressure characteristics of 
compacted silty soil. They have carried out 3 types of tests : 

a. Q tests : During these type of tests, compacted 
partially saturated soil samples are subjected to undrained 
triaxial shear with constant cell pressure. 

b. R tests : In these tests, compacted samples are placed 
in triaxial cell and are back saturated . Then, on these 
samples, undrained tests are carried out. 

c. S tests : These are slow drained tests on compacted 

back saturated soil. 

It may be noted that this silty soil has similar 
engineering characteristics as the present alluvium silt on 
which experiments have been carried out. 



4. Fredlund et al. (1978) have given the following 
equation : 

x = c'+(cr - u^) tan <p'+ (u - u ) tan <£ b (2.19) 

where <p D relates contribution of matrix suction to the 
shearing resistance. It may be noted that suction consists of 
two components — osmotic suction and matrix suction, (u - u^) 
is the capillary tension. Suction would be equal to the 
capillary tension (i.e. matrix suction ) only when components 
of osmotic forces were zero. This puts the limitation on degree 
of saturation range for which Eq. (2.18) is valid. 

5. Marsal (1979) suggested an approximate procedure for 
deriving the shear strength envelope for the partially 
saturated soil samples from UU test results. 

6. A more detailed analysis has shown that u^ would be 
very close to zero for high degrees of saturation (in general 
above 85% and in some soils above 90 % ) and that under these 
conditions the air would be occluded and the soil would be in a 
state called ' quasi-saturated' [Cruz et al. (1985)]. 

It is shown that effective stress envelope for partially 
saturated clay is not a straight line but consists of three 
straight line sections as shown in Fig. 2.5 depending upon the 
stress level. For first two sections, values of c'and <p' would 




be different for saturated and partially saturated soil while 
for third section, saturated and partially saturated test 
results would be aligned along the same straight line [Cruz 
et al. (1985) ] . 

7. Yudhbir and Korchoke (1987) have shown that within the 
narrow range of degree of saturation where air and water void 
continuity assumption needed for both Bishop et al. (1960) and 
Fredlund et al. (1978) approaches (Eqs. 2.18 and 2.19) is 
satisfied (for S range = 75 % to 85 % ) , both approaches are 
similar. They have also shown that Marsal's procedure to 
estimate <f>' from UU test results gives good results over wide 
range of degrees of saturation (even below 75 % degree of 
saturation where most of the field problems are encountered in 
residual soils ) . 

2.4 LITERATURE RELATED TO TRACK FORMATION DESIGN 

1. Rubin et al. (1970) determined the threshold 
stress-depth relationship for London clay by extrapolation from 
limited amount of data and superimposed it on maximum shear 
stress versus depth relationship as shown in Fig. 2.6 and the 
design criterion suggested was that for a subgrade of given 
threshold stress, the depth to the subgrade must be greater 
than that indicated by appropriate nominal axle load. 
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2. Agarwal and Yog (1975) proposed a new approach for 
design of substructure thickness based on effective stress 
analysis. Information about net state of stress was found out 
by adding up the stresses due to the live load and superimposed 
load (determined by Boussinesq's theory) and plot of total 
stress path of existing stress condition was developed. 
Assuming A = 0 (for compacted clay) and B = 1 for different k Q 
values, plots of effective stress paths of net state of stress 
were derived. For any soil, knowing shear strength parameter, 
k^' line can be drawn. Intersection of k^.' line with effective 
stress path (for net state of stress) would give the thickness 
of ballast and subballast layers. For cohesive subgrades, for 
different values of c', charts have been proposed between <p r 
and thickness of ballast and subballast layer. 

3. So (1978) carried out numerical analysis of track (Sec. 
2.2.2) and gave the family of charts. One family of charts for 
predicting maximum stress level induced in one component of 
track structure due to change in different track components of 
reference track was given. Thus, a series of families of charts 
was given for predicting the maximum stress level induced in 
different components of track structure. 

The family of design charts for subgrade stresses due to 
the vertical loading has been proposed (Figs. 3.65 to 3.67). 
Changes due to the tie-spacing, ballast depth and subgrade 
modulus of elasticity have been considered. Rail size and t^e 
support condition have insignificant effects on subgrade 



stresses. The design criterion being the induced stresses 
calculated from above design charts should not exceed the 
stress level required for fatigue failure determined by 
experiment. 

4. Tadatoshi et al . (A) have presented charts (Fig 2.7) 
from which depth of ballast and subballast can be determined in 
order to avoid the mud pumping. 

Mud pumping of roadbed is a phenomenon caused when roadbed 
soil is softened by the rainfall or groundwater ; the mud 
climbs up through the voids inside the ballast and pumped out 
on the ballast surface by a passing train load. Occurrence of 
mud pumping causes not only reduction of elasticity of ballast 
but also increases track subsidence, and as a result 
maintenance overburden is required. Occurrence of mud pumping 
is caused primarily by softening of roadbed soil by water, by 
use of unqualified soil, and by overloading to roadbed, 
according to field investigations and laboratory tests carried 
out by Japanese national railways [Tadatoshi et al. (A)]. 

From the empirical formula used by Japanese national 
railways (Sec. 2.2.1) , they presented the reduction of roadbed 
pressure achieved by change of track structure viz. narrowing 
the tie spacing [Fig. 2.8(a)], increase of ballast thickness 
[Fig. 2.8(b)]. They have also related upper roadbed pressure 
and mud pumping with liquid limit [Fig. 2.9(a)] and with CBR 
[Fig. 2.9(b)]. 
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Fig. 2.7 Decision diagram for roadbed depth 
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Fig. 2 . 8A Rel ationship between tie spacing end roadbed pressure 
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Fig. 2. oa Relationship between roadbed pressure, 
liquid limit, and mud pumping 
(Applied rail pressure is 5 t, and 
ordinary wooden tie is used.) 
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Fig. 2. OB Relationship between roadbed 
pressure, CBR, and mud pumping 
(Applied rail pressure is 5 tons and 
standard wooden tie is used.) 




5. Janin et al. (1983c) have argued that deflection of 
rail ( 2 sta ) calculated numerically will be lesser than the 
actual dynamic deflection because of presence of voids in 
actual dynamic condition. Now, as the vertical stiffness of 
rail remains same, i.e. 



( 2 . 20 ) 


Therefore, will be more than Q s ^. a * A procedure to 


calculate max Q^y n has been given from which max a ' ( ^ L y n 
found out by multiplying gt^ (obtained numerically) by 


can be 


f ” ax Q dyn 1 


where subscript 'dyn' indicates dynamic, and 

subscript 'sta' indicates static. 

In this way, curves in Fig. 2.10 have been generated. 
Then, an elaborate design procedure has been given which uses 
Fig. 2.10 and some other empirical relationships. The design 
criterion being the compatibility between dynamic stress on 
subgrade (Fig. 2.10) and the bearing capacity of the subgrade 
soil. 

This study observes that threshold stress (<x^ - cr^) design 
approach by Rubin et al. (1970) is based on fatigue behaviour 



- - - (o dyrf max) 

Maximum dynamic vertical stress exerted 



(e) Depth "ballast + sub-structure layer" 


Timber sleepers (length 2.60 m) 

Concrete sleepers (various lengths liom 2.20 m to 

2.60 m) 


FIG 2.10 Maximum dynamic vertical stresses exerted on the subq>-adf 



and guards against failure of subgrade but in practice, it is 
necessary that the foundation be stressed well below the limits 
of failure to prevent excessive track maintenance. 

6. Senac (1984) recognizes two important aspects in way of 
track design to resist vertical forces viz. track equipment and 
ballast formation considerations. 

Track equipment is defined as (I^/y^) / S where y^ is the 
distance between the neutral axis of rail and rail flange. 
Table 2.1 gives the track equipment standards adopted by SNCF 
(French Railways) for a formation of normal bearing capacity. 
The quality of formation has been classified based upon static 
rigidity, p (Table 2.2) and formation coefficient, C. 
Furthermore, SNCF specifies the following minimum substructure 
depth . 


30-55 cm under wooden sleepers. 

40-50 cm under concrete sleepers. 

depending upon soil characteristics, for tracks handling 
50000 t/day with high speed traffic and a minimum depth of 25 - 
3 5 cm on slower , less heavily trafficked line, subjected to 
favorable hydrogeological and substructure conditions. 


7 . Hall (A) has described the studies conducted for design 
of upgraded Richards bay coal export line in South Africa. 
Track structure parameters (rails, sleeper type and spacing) 



Table 2.1 Track Equipment standards (BY SNCF) 


Category 

Track equipment 

<Vy f >/5 

(TO ) 

Maximum 

axle load 

(t) 

Maximum Train 

speed 

(km/h) 

1 

100/140 

18 

70 

2 

141/179 

18 

85 

3 

180/235 

20 

95 

4 

236/350 

21 to 23 

115 to 160 

5 

> 350 

22 

200 


Table 2.2 Track formation quality (by SNCF) 


Formation Quality 

Static Rigidity 
(p in t/mm) 

Formation Coefficient 
(C in N/cm 3 ) 

Bad 

— 

20 

Poor 

0.6 

30 

Average 
or Acceptable 

1.7 

90 - 100 

Good 

atleast 3.5 

at least 180 












were specified beforehand. Hence, only formation design was 
required. Initially, a theoretical design based on mechanistic 
type analysis (often used for roadways) was evaluated. 
Simultaneously, a detailed study was carried out over a similar 
operating track. The theoretical design was then altered to 
give more practical final design based on the experience gained 
from the study of similar operating track. 

Failure criterion adopted for theoretical design was based 
on limiting the vertical strain. Depending upon track life 
required and traffic to be dealt, total number of repetitions 
in the track life (N r ) can be calculated. 

Then, a design track was assumed for analysis. Mechanistic 
model for traffic-induced fatigue cracking for stabilized 
formation was used as : 

N f = 10 9 * 1(1 " e/e b) (2.21) 

where = number of repetitions required for failure, 

e, = vertical strain at failure, and 
b 

e = maximum horizontal tensile strain at the base of 
stabilized layers (to be calculated from elastic 
or empirical models) . 

Thus, N f was calculated. For satisfactory design, N f > N r » 

Important findings of this study on existing track were as 


follows : 



a. Track defects show marked increase with the rainfall in 
the surrounding area. 

b. Mean moisture content vs. depth relationship remained 
same for all the tracks whether they were built in cutting or 
embankment provided water table is not intersecting the 
formation. This shows that moisture penetrates from top of 
track structure. The source of moisture being the rainfall. 

c. Mud pumping is independent of shear strength of 
formation and depends upon gradation of subgrade. Mud pumping 
occurs mostly on tracks built on clays and silts. 

2.5 LITERATURE RELATED TO OTHER IMPORTANT ASPECTS 

2.5.1 EFFECT OF TRAIN SPEED ON STRESSES 

During track analysis, speed effect is normally taken into 
consideration by a load augmentation factor based upon 
experiments. 

1. Austrian Method 

Austrian railways account for speed by a chart giving 
relation between load-increasing factor vs. speed [RDSO (B) ] . 

2 . German Method 

German Railways give the following formula [RDSO (B)]. 


O' 


speed 


cr t (l + s t ) 


( 2 . 22 ) 



where 


s depends on track and varies between 0.1 y to 0.3 y, 
y = 1 for speed 'V ' s 60 kmph, 

= 1 + - - f40 60) ■ for V > 60 kmph, and 

t = 3 for a statistical safety of 99.7 %. 


3. Japanese Method 

According to Japanese practice, speed impact factor can be 
found out from empirical relation [RDSO (1987)] given below : 


i = 1 + 0.5 V / 100 (with joint rail) 
i = 1 + 0.3 V / 100 


(2.23) 


'i ' should not exceed 1.8 


2.5.2 CHARACTERIZATION OF TRACK FOUNDATION 




km 
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There are several ways in which elastic foundation can be 
characterized as given below. 


1. Formation coefficient 'C' (kg/ cm 3 ) [RDSO (B) and Senjac 
(1984) ] . 

2. Modulus of elasticity of rail support or modulus of 
rail-road reaction 'k' (kg/cm 2 ) [Talbot et al. (1923)]. 

3. Static Rigidity 'p' (kg/cm) [Senac (1984)). 


The modulus of elasticity of rail support 'k' had been 
introduced as a measure of vertical stiffness of rail support 
[Talbot et al. (1923)]. It may be defined as the pressure per 



unit length of each rail required to deflect the track per 
unit. 


k 


equivalent load / tie 


Average depression x spacing 


It represents the stiffness and yieldability of tie, 
ballast and roadway but does not involve the stiffness of rail. 
It is related to tie-spacing, tie-dimensions, depth of ballast, 
quality of tie, solidity of roadway, character of tamping and 
surface and other conditions of the track and its maintenance. 
Naturally, this value must be determined by tests at different 
locations and over a wide range of conditions. It is denoted by 
'u' or 'k'. 

In Z immermann ' s method, either ' C or 'k' is used and it 
is assumed that this 'C' or 'k' characterizes the continuous 
foundation [RDSO (B) ] . This assumption does not have much 
effect and Z immermann' s method is good to use because although 
the cross-ties constitute a discontinuous support, their width 
forms a considerable part of the rail length, and , as the 
distance from a single wheel load to the point where the track 
depression is zero covers from 4 to 6 tie spacings, the 
assumption of continuous support may not be expected to be 
generally much in error for usual conditions of good track, 
especially as the method used for determining the value of 
modulus of elasticity of rail support 'k' will carry into this 



value some of the conditions attending the discontinuous 

support . 

Accuracy of analysis using 'k' may be affected because / k / 
may not be a constant because initially, some load will cause 
readjustment of pores and gaps in the track structure. 
Conditions of track may be classified into following 3 groups 
as shown in Fig. 2.11 [Yudhbir et al. (1991)]. 

a. Condition 'A' represents 'Well tamped' track. 

Load-settlement curve for this track can be represented by 

straight line passing through the origin. 

b. Condition 'B' represents 'track condition where bending 

of tie is necessary before full bearing is obtained'. 
Load-settlement curve for this track can be represented by 

straight line which does not pass through the 

origin. 

c. Condition 'C' represents 'poor' track. Load— settlement 
curve is not a straight line. 

Concept of u^ and u^ : 

The load method for the calculation of track modulus 
[Kapoor et al. (1967a)] is given below 

(10.7 Q.) 4/ ^ 3 

Initial range u. = yy-r ~rr\ lbs/ in/ in (2.24) 

1 <I r 7 * y sl> 7 




Fig. 2.11 TRACK DEPRESSION VS. LOAD RELATIONSHIP 


FDR DIFFERENT TRACKS. 



where 


Q = total vertical load 


Q 1 = initial load required to bring the track 
in elastic range. 

y sl = deflection of sleeper under load Q. 
u^ = modulus of railroad reaction in initial 
range. 


elastic range u g = 


(10.7) 4//3 (Q - Q x ) 4/3 

<y s -y s i)3 4/J 


lbs/in/in (2.25) 


Note 1 kg/ cur =6.89 psi 


where 


u = modulus of railroad reaction in elastic 
e 

range. 


Bearing area of sleeper [Mundrey (1988)] : 

On empirical basis, the bearing area below sleeper which 
takes the load is as shown in Fig* 2.12. It is assumed that 
stresses are uniformly distributed on bearing area. 

Bearing length / sleeper as given by Mundrey (1988) — 

131.42 cm (as shown in Fig. 2.12). 

Bearing length / sleeper as given by Saxena and Arora 

(1987) = 180 cm. 

Bearing length / sleeper as calculated from data given by 
Senac (1984) = 166.2 cm (adopted here). 
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Fig. £,1£ BEARING AREA DF SLEEPER 




Relation between 'C', 'p' and 'k' [see Yudhbir et al. (1991)] : 


k = b c (2.26) 

Bearing Area per sleeper ' F' 

where b = 3 

(2 x spacing 'S') 

C = p / F , and (2.28) 

P = R / y t ( 2 . 29 ) 

with R = Reaction of sleeper, and 

y^. = Vertical deflection of sleeper due to force 'R' 

p = 2 k S (2.30) 

2.5.3 SOIL IMPROVEMENT METHODS 

This can be achieved by methods given below [RDSO (C) and 
Tadatoshi et al. (A)]. 

1. Soil stabilization : A list of usual methods for 

roadbed stabilization has been given in Table 1.1. 

2. Soil drainage : Proper drainage plays very important 

role in controlling the mud pumping. 

3 . Soil compaction : Proper compaction increases the 
threshold stress for soil. 

4. Soil reinforcement : Japanese railways have suggested 
that for high speeds, heavy axle load and increased traffic 



density, soil reinforcement may be used in addition to the 
usual thickness of ballast as obtained from theory of 
elasticity [RDSO (C) ] . 

2.5.4 FIELD MEASUREMENTS AND MODEL TESTING 

1. Field measurements by Indian railways [RDSO (A)]. 

Empirical formulae need the information on track 
modulus and hence to find this out for Indian tracks, field 
studies were carried out in 19 60 ' s at Ambala - Ludhiana section 
which resulted in civil engg. report no. C 47, Nov 1966. Later, 
track modulus (initial and elastic) values had been recommended 
as given in civil engg. report no. C-67, Sept. 1968 of RDSO 
(Research, Design and Specifications Organization, Indian 
railways) . 

Track with wooden sleepers : u^ = 50 kg/cm/cm 

u g = 250 kg/cm/cm 

Track with metal sleepers : u^ = 75 kg/ cm/ cm 

u e = 300 kg/cm/ cm 

Till 1960, no field tests were conducted to check whether 
empirical formulae apply to Indian track conditions. Hence in 
1962-66, dynamic field tests were conducted to measure the 
maximum pressure on top of formation which resulted into civil 
engg. report no. C-47, C-48, C-55 and C-56. 

Since dynamic field tests were conducted only on tracks 
with N+3 sleeper density and the track structure later 



consisted mainly of N+6 sleeper density, and dynamic tests 
posed problems due to line capacity, static tests were 
conducted by RDSO in 1970-74 to find out pressure distribution 
along the depth. The results were presented in civil 
engg. report no. C-114, and C-131 of RDSO. 

Long haul trains having 4500 t and 9000 t trailing loads 
were introduced in 1980's on Indian railway tracks. To 
ascertain the effect of heavy haul trains on the stress coming 
on the formation, dynamic tests were carried out by RDSO. 
Finite element analysis requires elastic modulus of formation 
soil as one of the input parameters. Hence, RDSO in 
consultation with IIT Kanpur conducted cyclic plate load test 
and wave propagation test at Moradabad division, northern 
railway [Kumar et al. (1990), Kumar et al. (1991a) and Kumar et 
al. (1991b)]. For these tests, in general, test procedure 
described in IS 1888-1982 and IS 5249-1977 were followed. 

2. MODEL TESTING BY ORE 

Janin et al. (1982) undertook full scale laboratory 
testing at Vienna arsenal to determine the settlements and 
stresses in the track structure under static and dynamic 
loading. These tests were conducted in a large rectangular 
reinforced concrete box having a section of the track 
foundation structure. This substructure was loaded using a 
single sleeper with 90 kN load. Displacements and stresses were 
measured at the interface of all layers. 



Tests were conducted over four model track beds having 
varying thickness of substructure layers. Subgrade was made up 
of clay of low compressibility. Each test was carried out for 6 
million cycles with application of approximately sinusoidal 
loading waveform from 0 to 90 kN at 10 Hz frequency. 
Simultaneously, computer program 'ROSALIE' was used to 
calculate displacements and stresses for these tracks. Measured 
displacements and stresses at first loading for track beds have 
been compared with calculated results from 'ROSALIE' . 

At Vienna Arsenal, Janin et al. (1983a) further carried 
out similar box tests on a different track bed (having no 
gravel layer below ballast) . Subgrade was made up of sandy silt 
of low compressibility. 

However, most detailed testing programme was carried out 
at Derby [Janin et al. (1983a)]. An area of 20 m x 4.5 m 
containing track panel was selected for test site. One half of 
this site was filled with soft boulder clay and other half with 
stiff keuper marl clay as subgrade material. Instrumentation of 
this site consisted of isolating the area from surrounding by a 
polythene membrane, pressure cells, accelerometers, 
potentiometers and LVDT recordings of displacement response at 
the interface of different layers and hydraulic piezometers 
just below the clay surface for pore water pressure 
measurement. Dynamic and static tests consisted of loading rail 
head and single sleeper by servo controlled hydraulic vibrator. 
Rolling tests were carried out by propelling a two axle 
vehicle (16t axle load) across the test site. 



Four different depths of formation were investigated at 
Derby. Measurements were taken during rolling, static and 
dynamic loading to evaluate the effect of depth of formation on 
the following parameters. 

1. Elastic and permanent settlement of ballast, 
subballast, blanket and subgrade. 

2. Vertical, longitudinal and transverse pressures at the 
subgrade surface. 

3. Pore water pressures in the subgrade. 

2.6 SUMMARY 

The present study to predict the behaviour of soils under 
the railway formations subjected to cyclic loading requires 
determination of induced stresses on the subgrade by 
appropriate mathematical model and the evaluation of threshold 
stress of subgrade soil by geotechnical testing. 

It has been observed that before 1980, induced stresses 
were determined by empirical formulae. With high speed traffic, 
sophisticated mathematical computer oriented models were 
increasingly used. These models can be broadly classified under 
the finite element models and models which use Burmister's 
solution for multilayer system [Burmister (1945)]. These all 
models give the solution by assuming equivalent static load 
which is obtained by augmenting actual wheel load by speed 


factor. 



It is observed that there are very few simple and yet 
accurate models available to predict induced stresses. Also, it 
is seen that most of the two dimensional models predict higher 
stresses as compared to the measured stresses. In the present 
study, a simple model 'ELASTIC. F' has been developed and it is 
shown that this model gives quite accurate results. Also, two 
dimensional finite element model '2D8N' has been presented 

along with a modification which enables this model to predict 
quite reliable vertical stresses. In addition to this, a three 
dimensional linear elastic model '3D20N' has been developed 

which predicts result quite similar to 'MULTA' model which is 
like a closed form solution. 

In this literature review, it is found that , in general, 
there are two different approaches to evaluate the threshold 
stress for soils subjected to cyclic loading. Following Sangrey 
(1968) , critical level of repeated stress (similar to threshold 
stress) has been evaluated from stress path movement during 

undrained cyclic loading and is taken as stress level 

corresponding to end of equilibrium line. Also, Rubin et al. 
(1970) have observed that plastic strain development with 
cyclic loading shows terminating and non-terminating behaviour 
for undisturbed samples depending upon cyclic stress and based 
on this, threshold stress can be evaluated. Guidelines given by 
Sangrey et al. (1978) to evaluate the threshold stress requires 
determination of k . and hence not easy to use. In the present 
study, it has been shown that threshold stress for compacted 
quasi saturated soils during undrained triaxial tests can be 



determined by both above mentioned approaches and they yield 
same results. Also, a generalized relationship between and 
threshold stress ratio has been proposed. 

Literature review also reveals that relationship between 
cumulative plastic strain and number of loading cycles is 
available for undisturbed London clay [Rubin et al. (1970)]. 
However, this model requires determination of elastic strain, 
e e (measure of cyclic stress ratio) and hence cannot be used 
directly. In the present study, a generalized model for plastic 
strain development with number of loading cycles has been 
proposed. 

In addition, it may be noted that data available on 
compacted samples subjected to cyclic loading is rather scarce. 
The present study presents and analyzes the data (both static 
and cyclic) on compacted samples. Also, an attempt has been 
made to explain the static data on compacted samples in the 
critical state framework. 

Literature review also reveals several design methodology 
used by different railways in the world. However, they cannot 
be directly used in Indian conditions. In the present study, a 
methodology for design of railway formation has been 
proposed. 



results obtained by other models and field studies is given in 
Sec. 3.5. Section 3.6 compares behaviour of these three models 
in terms of accuracy, results obtained and CPU time. 

3.2 RIGOROUS ELASTIC METHOD (ELASTIC. F MODEL) 

3.2.1 INTRODUCTION 

This model consists of determining the stresses at the 
sleeper ballast interface by considering the rail as a beam on 
elastic foundation. The stresses inside the track structure are 
then determined by converting layered system into equivalent 
homogeneous system by Odemark method and using Newmark's 
charts. This method has been described below in detail. The 
parametric studies have been carried out using this method. 


3.2.2 THEORETICAL BACKGROUND 

The governing equation for beam on an elastic foundation 
is given as 


^bm ^bm 


d I * * 4 y 


-k y + g (X) 


(3.1) 


where 


I = Young's modulus and moment of inertia of 

'bm ' bm 

beam, 

x = distance of point where deflection is 
required from point where load is acting, 



y = deflection of beam, and 
q(x) = externally applied load. 

For the simplest case of single concentrated load / Q / 
acting in the vertical plane of symmetry on an infinitely long 
rail, Timoshenko and Langer (1932) have given following 
solution : 

y r = Qt 7(£)/2 k 1^ (3.2) 

(3.3) 

(3.4) 

(3.5) 

P s is the stress below the sleeper where # y ' has 
been calculated. p g is assumed to be uniformly distributed over 
the bearing area per rail seat (Sec. 2.5.2) of sleeper referred 
here as Load area. 

and 7}(0 = e"^[cos (?) + sin (£) ] (3.6) 

[Hetenyi (1961)] 

In Eq. (3.6), 


where 

y r = Q X [7 J(e i ) + U(? 2 ) ] / 2 k Lh 




5 = x / Lb 


(3.7) 



Referring to Fig. 3.1, for £ 1# x in Eq. (3.7) will be the 
distance from load Q 1 and for £ , x will the distance from load 

V 

This method to determine stress at sleeper ballast 
interface is also known as Zimmermann's method [RDSO (B) ] . 

Track structure is generally made up of 3 layers viz. 
ballast, subballast and subgrade soil (Fig. 1.1). Odemark 
method is extensively used for vertical stress calculation in 
multi-layered systems [RDSO (B) ] . By Odemark method, layered 
system can be converted into homogeneous one (Sec. 2.2.1). 
Odemark method gives good results only when soft layer is 
underlain by stiff layer. 

Referring to Fig. 2.1, the homogeneous equivalent depths for 
ballast and subballast layers are given as 

d b - °- 9 d b /v*I 

d sb = 0-9 d sb ~'J E sfc/ E s 

Newmark's charts [Jumikis (1966)] can now be applied to 

: k 

find out the stresses induced at this new equivalent depth (d fc + 
d*^) from the stresses at sleeper ballast interface. 

3.2.3 DESCRIPTION OF RIGOROUS ELASTIC METHOD (ELASTIC. F) 

Sleeper and Load configuration for double wheel load case 
are shown in Fig. 3.2. Load tie is a tie directly above which, 


(3.8) 

(3.9) 
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Fig. 3.2 SLEEPER AND LOAD CONFIGURATION 
FOR DOUBLE WHEEL LOAD CASE 



the wheel load is acting. In this method, load tie and 4 ties 
on either side of it, are considered for stress calculation. 


Step 1 - Calculate the stresses at sleeper-ballast interface as 
given below : 

If track characterization is given in terms of formation 
coefficient, C or static rigidity, 'p' , convert it into 'k' as 
given below (Sec. 2.5.2) : 

P 

k = (3.10) 

2 S 

and k = b .C (3.11) 

Uniform stress on bearing area per rail seat of sleeper 
(load area) can be obtained by Eqs. (3.4) and (3.5). 

Thus, stress in any load area due to load at any place can 
be obtained. 

In this analysis, interaction of four sleepers on either 
side of load sleeper has been considered (Fig. 3.3). Due to 
symmetry, only one quarter of track need to be analyzed (Fig. 

3.4). The stresses on load areas LI to L5 are taken equal to 

summation of stresses due to loads Q 1 and Q 2 while stresses on 
load areas L2 ' and L5' are taken equal to stresses due to load 
only. This is because it has been observed [ (Talbot et 
al. (1923) ] that effect of load Q2 need to be taken only on 

four sleepers on left (i.e. on LI to L4) and on sleeper L5 for 

calculating stresses below point A (Fig. 3.4). Stresses on load 
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Fig. 3,3 Assumed load distribution below the. sleeper 
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Fig. 3,4 Assumed load distribution below the sleeper 


areas L6 to L10 and L7 ' to L10' can be obtained by symmetry. 

Step 2 - Conversion of layered track structure into homogeneous 

system by Odemark method (Fig. 2.1) : 

. * * 

Equivalent heights d fc and d^ are calculated from Egs. (3.8) 
and (3.9), 

Step 3 - Calculation of stresses in track structure : 

Referring to Fig. 3.3, stresses at any point in the track 
structure can be found by superimposing the effect due to load 
areas LI to L10 and L2' to L10'. 

Stresses due to any rectangular load area (length a, width 
b) at any point in the track structure can be found by 
Newmark's solution [Jumikis (1966)] as given below : 

Vertical stress at corner of any rectangular area due to 
uniformly distributed load on total area is given by 

2 in n ( m 2 + n 2 + 1) 1//2 in 2 + ii 2 + 2 
— x — 

m 2 + n 2 + in 2 n 2 + 1 in 2 + n 2 + 1 

2 in n ( m 2 + n 2 + 1) 1//2 

+ tan — — 

in 2 + n 2 - in 2 n 2 + 1 


4 71 


(3.12) 



where 


(3.13) 

(3.14) 


with 


® = a / z d 

s - b / Z(3 

z 3 = depth of point from surface at which stresses are 
required. 


or, 


( cr z ) c = q 


n 


(3.15) 


where 

q^ = intensity of pressure on rectangular load area. 
k n = influence factor to be taken from charts given by 
Newmark. 


Eqs. (3.12) or (3.15) can be used for calculation of 
vertical stress at a point situated either inside or outside 
the rectangular area by the method of superposition [Jumikis 
(1962) ] . 

A computer program based on this method has been written 
[Yudhbir, Kameswara Rao and Shahu (1991)]. 


3.2.4 LIMITATIONS 

1. This model requires experimental determination of 

modulus of railroad reaction, k or formation coefficient, C in 

addition to E , E. , E . and E . 

r' b' sb s 

Although, no relationship between E and the value of 

b 



modulus of railroad reaction (k) or formation coefficient 'C' 
is available in literature, modulus of subgrade soil does have 
bearing upon the value of modulus of railroad reaction (Sec. 
2.5.2). For a rough estimate, the following procedure may be 
adopted : 

Senac (1984) gives the track formation quality in terms of 
p and C (Table 2.2). Bowles (1984) suggests a range of Young's 
modulus for different soil types (Table 3.1). Table 3.2 shows 
approximate correlation between C and E g based on the 
assumption that formation quality of track is determined by the 
type of subgrade on which it is laid. 

2 . Odemark method should not be used for systems where 
stiff layer is underlain by soft layer. This limitation also 
applies to this model. 

3 . This model can not simulate factors such as tie 
characteristics, outer track geometry, boundary conditions etc. 

3.2.5 COMPARISON OF RESULTS 

Comparison of results has been given in Sec. 3.5. 

3.2.6 RESULTS AND DISCUSSION 

TRACK 1 has been considered as reference track. TRACK 1 is 
the representation of typical Indian track. Details of this 
track are given in Appendix Al. Parametric studies have been 



Tabic 3.1 Range of Young's modulus values for different soils 
[Bowles (1984)] 


Soil type 

E s 

(MPa) 

Very soft 

2-15 

Soft 

5-25 

Medium 

15 - 25 

Hard 

50 - 100 


Table 3.2 Empirical Correlation between formation coefficient, C 
and Young's modulus of subgrade, E 

S 


Formation quality 

C 

Adopted for parametric studies 

c 

E s 

(N/cm 3 ) 

(N/cm 3 ) 

(MPa) 

Bad 

>20-30 

25 

10 

Poor 

>30-90 

75 

20 

Average 

> 90 - 180 

125 

30 

Good 

> 180 

180 

50 






carried out by varying range of parameters with reference to 
this track. Range of parameters used are as follows : 

Depth of subballast = 0.15, 0.35, 0.55, 0.75 m 

Young's modulus of subgrade soil = 10, 20, 30, 50 MPa 

Moment of inertia of rail = 1600, 2158, 3055 cm 4 

Tie spacing = 0.50, 0.65, 0.76 m 

The results of parametric studies are given in Table 3.3 

Figure 3 . 5 shows the variation of vertical pressure on top 

of subgrade under the rail seat (cr with Young's modulus of 

subgrade (E ) . It will be shown later (Sec. 3.3.7 and Fig. 
s 

3.39) that the maximum vertical pressure on top of subgrade 
occurs below the rail seat and thus, cr ^ represents maximum 
vertical pressure on top of subgrade. 

From Fig. 3.5, it may be seen that as the soil becomes 
stiff, <r increases. Also, as the total depth of formation 
(d^.) decreases, cr ^ increases. It can be seen that effect of 
soil stiffness on cr ^ is more when total depth of ballast and 
subballast is less as expected. Figure 3.6 shows the variation 
of or with total depth of formation for four soil types. The 
trend in Figs. 3.5 and 3.6 is same as obtained by So (1978). 

Figures 3.7 and 3.8 show the variation of vertical stress, 
cr under the wheel load inside the subgrade soil up to 1 meter 
depth for varying total depth of formation for soft soil (E g = 
10 MPa) and stiff soil (E = 50 MPa) respectively, cr decreases 
by only 30 % (Fig. 3.7) for d t = 1.0 m while for d t = 


0.4 m , 
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75 



49. 56 
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File 

Name 

E b 

CMPa3 

E sb 

CMPaD 

E 

s 

C MPaD 

d b 

C crrO 

d sb 

CcnD 

I 

r 

Ccm ) 

Tie spacing 

CcnO 

c on top of subgrade 

CkPaD 

B2T1 

150 

80 

20 

25 

35 

21 58 

50 

54 74 

B2T3 
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01 
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50 

25 

15 

2158 
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D1SB3 


100 
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DSSBls 
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00 
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2158 
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100 

50 

00 

60 

2158 

65 
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Fig. 3.7 Variation of CJv inside the subgrade 
for E s =10 MPa for different df 
(ELASTIC- F model) 



82 


| CT V (kPa) 



Fig. 3.8 Variation of OV with depth for different 
dt (ELASTIC-F model) 




the corresponding decrease is 50 % (E = 10 MPa) . Thus, it can 
be seen that from point of view of stress dispersion, subgrade 
soil is more important if depth of formation is less. 

Figures 3.9 and 3 . 10 show the variation of & with depth 
under the wheel load for different soil types for d t equal to 
0.4 m and 1.0 m respectively. 

From Figs. 3.7 to 3.10 it can be seen that for a given 
soil type, provision of more total depth of formation than a 
certain limit does not reduce the vertical stress on top of 
subgrade under the rail seat by significant amount. 

Figure 3.11 shows the variation of maximum vertical stress 
on top of subgrade for different proportions of ballast and 
subballast depths but keeping the total depth of formation 
equal to 1.0 m. 

From Fig. 3.11 , it is clear that increase in depth of 

ballast layer but keeping total depth of ballast and subballast 
layer constant does not reduce cr^ significantly. 

It is to be noted that while using ELASTIC. F model, 
empirical correlations of ' C' and 'E ' values for similar types 

b 

of formation and soil have been used (Sec. 3.2.4) . Hence, it 

is useful to check the variation of <r ^ values with formation 
coefficient ' C '. For d. = 0.4 m and E = 10 MPa, 'C' was 

u S 

3 3 

changed from 25 MN/m (recommended value) to 250 MN/m (extreme 

value in all categories) . The change in maximum vertical stress 

on top of subgrade is only from 55 KPa to 66 KPa (Fig. 3.12) . 

3 3 

Next, 'C' was changed from 25 MN/m to 75 MN/m (recommended 
value for second category) . The change in maximum vertical 
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cr v (kPo) 



Fig. 3.9 Variation of <T V with depth for different 
E s (ELASTiC-F mode! ) 
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f OV (kPa) 



Fig. 3.10 Variation of CTy inside the subgrade 
for dt =1.0m for different E s 
(ELASTIC-F model) 




OVtCkPa) 3! CT vt (kPa) 



db (m) 


.3.11 Variation of CT v t for different proportioning 
of ballast and subballast depth for dt = 1.0m 
(ELA5TIC-F model) 



C (N/cm 3 ) 

Fig. 3.12 Variation of CTvt with formation coefficient 1 ^ 
for dt =0.4m (ELASTIC*F model) 




stress on top of subgrade is only of the order of 4 KPa (7 %) . 

Also it will be shown in Fig. 3.41 that 8 (or C) vary 

significantly only with E s while other track parameters and 

depth of formation has negligible effect on 5 (or C) . Thus, the 

correlation between ' C' and 'E ' (Table 3.1) is fairly good. 

s 

Figure 3.13 shows the effect of tie spacing on maximum 

vertical stress on top of subgrade. As the tie spacing 

increases, <r ^ increases as also noted by Stewart and Selig 

(1982). Two extreme values of tie-spacing have been used (0.50 

m and 0.76 m). It may be seen that tie spacing has small effect 

on <r vt for higher d^. (for d^. = 0.60 m, effect is 5 % while for 

d t = 0.40 m, it is 15 % ), in general , but this effect is even 

less for soft soil as compared to stiff soil (negligible for E 

= 20 KPa ; 6 % for E = 50 KPa ) . 

s 

Figure 3.14 shows the effect of moment of inertia of rail 
on vertical stress on top of subgrade under the rail seat for 
depth equal to 0.4 0 m. It may be noted that rail type has 
insignificant effect on cr 

Figures 3.15 and 3.16 show the effect of Young's modulus 
of ballast (E fc ) and subballast (E gb ) on cr vt respectively. The 
Young's modulus values for both ballast and subballast layers 
within their given range have insignificant effect on cr^. 
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, 3.13 Variation of Ovt with tie spacing for 

different track conditions (ELASTIC *F 
model ) 



Fig. 3.14 Variation of <7 v t with moment of inertia 
of rail , lr 
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Fig. 3.15 Variation ot (Xvt with Young's modulus 
ot ballast (ELASTIC F model ) 
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Fig. 3.16 Variation Ot dyt with Young's modulus 
of subballast ( ELASTIC -F model) 





3.3 TWO DIMENSIONAL FINITE ELEMENT ANALYSIS 


3.3.1 INTRODUCTION 

Finite element method is a well established technique and 
has been described in a number of books [e.g. Zienkiewicz and 
Taylor (1989); Cook et al. (1989); Rao (1989); Reddy (1986)]. 
Hence only the highlights of the routines developed have been 
described here. 

All the finite element routines are based on displacement 
formulation of elasticity equations. All routines use linear, 
elastic material properties. Two dimensional finite element 
analysis of track has been carried out in order to calculate 
deflections and stresses inside the subgrade. Two dimensional 
finite element model for railway track has been described below 
in detail. Stress distributions and displacements for a typical 
track have been presented and a parametric study was carried 
out using two dimensional finite element analysis with 
reference to this track. 

3.3.2 DESCRIPTION OF TWO DIMENSIONAL FINITE ELEMENT ROUTINES 

Two finite element routines are developed , each of which 
has capabilities to solve both plane stress and plane strain 
problems. Names given to these routines are : 



1. 2D4N 


2. 2D8N 

2D4N is a finite element routine which uses 4-noded linear 
isoparametric rectangular elements. 2D8N is a finite element 
routine which uses 8-noded parabolic isoparametric rectangular 
elements. Both these programs use similar formulations, 
techniques, inputs and outputs. Vertical and lateral 
deflections and stresses, and shear stresses at all nodes can 
be calculated. 2D4N routine is equipped with semi-automatic 
mesh generation facility while 2D8N routine is equipped with 
both semi-automatic and automatic mesh generation modules. 

Both routines were tested for standard cantilever beam 
problem [Fenner (1986)] and Boussinesq's line load problem 
[Harr (1966)]. It was observed that 2D8N module gives better 
results with fewer number of elements as expected and hence, 
further analysis was done with 2D8N only for two dimensional 
finite element analysis, comparison for Boussinesq's line load 
problem has been given in Appendix Bl. 

2D8N routine can handle layering of different materials 
and all types of boundary conditions. The routine has facility 
to simulate joint between different layers by 6— noded linear 
joint elements of zero thickness [(Goodman et al. (1968); 
Burgohain and Shah (1977)]. The results obtained from this 
routine have been checked with the standard matrix given by 
Vijayan (1981) for interface element . 

This routine has the facility to simulate elastic 



continuum by static infinite elements [Curnier (1983)]. 
However, use of infinite elements needs special kind of mesh 
discretization. Mesh should spread out in radial direction from 
the load element [Zienkiewicz et al. (1983)]. It has been noted 
that infinite elements with rectangular mesh discretization do 
not give better results. As the track problems do not allow 
this special kind of radial mesh discretization due to the 
presence of layering, infinite elements were not used for track 
problems. 

Input in the automatic mesh generation module consists of 
basic parameters (e.g. total number of nodal points, material 
types, elements etc.), material properties, total number of 
compartments in x and z direction and length of each 
compartment, boundary conditions at four faces etc. Mesh 
generated will have rectangular lay out. 

For generation of meshes of any shape, semiautomatic mesh 
generation facility can be invoked which requires input as 
basic parameters, material properties, nodal point data, 
element connectivity data, boundary conditions, etc. 

Output from the routine consists of both horizontal and 
vertical stresses, shear stresses and principal stresses at all 
Gaussian points and extrapolated stresses at all nodal points 
along with displacements at all nodal points. Stress 
extrapolation is carried out using linear shape function [Cook 
et al. (1989)]. However, to cope with the enormous amount of 
data, a post processor facility has been provided using which 
important output can be obtained in the required form. Thus, 



for e.g. , output related to deflection and stresses at nodal 
points along depth below the load or in any horizontal 
compartment at given depth can be obtained. 

As already mentioned, routines were checked for plane 
strain Boussinesq's problem with line load [see Appendix B1 ; 
Yudhbir , Kameswara Rao and Shahu (1991)]. An interesting 

feature was noted at this time. For mesh with 250 elements, and | 

{ 

bottom and lateral boundaries at 6.0 m from the load, it was 
observed that computed vertical stresses and displacements were 
as good as exact but lateral stresses did not match with exact 
results throughout length and depth of mesh. Moreover, lateral 
stresses were found to vary with Poisson's ratio. Mesh 
refinement and change in boundary conditions were not much 
helpful. However, it was found that when lateral and bottom 
boundaries were fixed at 130 m, lateral stresses match with j 
exact results and do not vary with Poisson's ratio anymore. 

The notations used in program are shown in Fig. 3.17. 

3.3.3 TWO DIMENSIONAL FINITE ELEMENT MODEL FOR RAILWAY TRACK 
(2D8N MODEL) 

Reference axes chosen for finite element models are shown 
in Fig. 3.18. 

Initially, both plane stress and plane strain cases were 
tried. For plane stress case, longitudinal section of track was 
analyzed in the similar way as for Illitrack model [Selig et 
al. (1979)]. Longitudinal section of track is a section which 
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Fig, 3,18 TRACK STRUCTURE AND THE REFERRED DIRECTIONS 
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runs along the rails while transverse section is one which runs 
across the rails i.e. along the sleeper (Fig. 3.19). 

Stresses calculated from plane stress case for a 
particular point at a given depth in longitudinal section, 
would be the average uniform stress at that depth in the 
transverse section. For example, if stress at the point A (at 
depth z) was found out as 6.5 KPa by plane stress analysis then 
uniform stress distribution at depth z in transverse direction 
b e 6.5 KPa. But it is well known that stresses in 
transverse section at any depth are non uniform [for e.g. , see 
Selig et al. (1979)]. Also, the thickness of longitudinal 
section is varying from top to bottom and there is no rational 
way to decide which thickness is to be taken for analysis and 
if the thickness is to be varied from top to bottom then at 
what rate with the depth it should vary. Hence the plane stress 
case was subsequently discarded. 

For plane strain case, transverse section of railway track 
is considered. For this case, all the sections in y-direction 
(Fig. 3.18) must be the same in loading and geometry everywhere 
for a very long distance so that strain in y-direction is zero 
i.e. 0 ^= o and need not enter the analysis. For railway track, 
there are two difficulties in meeting with this requirement. 
First, the sleepers are not continuous but are discrete in 
y-direction. Second, due to the presence of concentrated wheel 
load, loading symmetry is not fulfilled. 

Though sleepers are discrete in y-direction, they are 
assumed to be continuous in that direction but the overall 
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Fig, 3,19 TRACK SECTIONS 




flexural stiffness El of the continuous sleepers was made equal 
to that of the discrete sleepers. 

For example - Considering the sleeper section as shown in 
Fig. 3.20, moment of inertia of discrete sleeper is given as 

I tl = b t3 / 12 (3.16) 


where ^tl ~ moment of inertia of discrete sleeper. 

Keeping the same sleeper thickness / t / and width 'b' for an 
equivalent continuous system, 


E tl I tl E t2 I t2 

(3.17) 

b t 3 E tl / 12 = St 3 E t2 / 12 

(3 . 18) 

r .. ' 

E t2 " b E tl/ E 

(3 . 19) 

where 


= modulus of elasticity 

E t2 ~ moc ^ ulus of elasticity 
I t2 = momen t of inertia for 

for discrete sleeper, 

for continuous sleeper, and 

continuous sleeper. 

In order to satisfy loading symmetry for plane strain 

case, wheel load is converted into equivalent line load. The 

magnitude of the line load can be decided by the following two 


empirical ways : 


1. Figure 3.21 shows the wheel configuration for a typical 
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Fig. 3.20 SLEEPER SECTION 
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Fig, 3.21 WHEEL CONFIGURATION 






train [Selig et al. (1979)]. The spacing between different 
wheels of rail car is shown there. To be on the safer side, 
taking the spacing between two wheel loads equal to 1.70 m for 
the purpose of converting the wheel load into line load, 
average line load intensity can be obtained as : 

For TRACK 2, load of one wheel = 144.657 kN 
Then, equivalent line load on sleeper = 144.657/1.7 kN/m 

= 85.092 kN/m 

This method gives the average line load. 

2. Load distribution given by Talbot et al. (1923) for 
wheel placed in the middle of two sleepers is shown in Fig. 
3.22 and the load distribution given by Japanese railways for 
wheel on top of sleeper [Tadatoshi et al. (A)] is shown in Fig. 
3.23. 

Superposition of the above two cases gives required load 
distribution as shown in Fig. 3.24. Maximum line load intensity 
can be calculated as below (Fig. 3.24) : 

Maximum load is taken by sleeper beneath the load (i.e. 
load sleeper) = 0.43 Q 
Width of sleeper = b 

= 0.254 m 

Hence, line load intensity for load sleeper 

= 0.43 x 144.657/ b 


= 244.9 kN/m 
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Fig. 3.24 LOAD DISTRIBUTION — PRESENT PROBLEM 



Thus, maximum line load intensity 


= 244.9 kN/m 


Average line load intensity can be calculated as follows : 


Considering range 1 [(the major load areas), Fig. 3.24] : 


(0.43+0. 27+0. 25+0. 25+0. 25) Q 
Average line load intensity = 


5 x 0.495 


= 84.75 kN/m 


Another difficulty in modelling occurs due to the presence 
of rail section. It is due to many reasons. Firstly, for two 
dimensional analysis, sleepers are assumed to be continuous in 
longitudinal direction. Thus, if rail element is modelled as 
two dimensional (2-D) linear element, flexural action of rail 
will be prevented and load distribution will be through direct 
shear action. This will induce error in results. 

Secondly, while actual contact area between discrete 
sleeper and curvilinear rail is known, contact area between 
assumed continuous sleeper and rectangular rail section will 
have to be assumed without any rational basis. 

Thirdly, the determination of magnitude of equivalent line 
load in place of concentrated wheel load on top of rail head 
will be very difficult. 
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Fig, 3,25 Track formation originally used for calculations 
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Fig. 3.26 Discretization for mesh A (2D8N model) 
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Fig. 3.27 Discretization for mesh B (2D6N model) 
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Fig. 3.28 Discretization for mesh C (2D8N model) 


mesh B and mesh C , data for TRACK 2 [ (same as FAST track 
section, Selig et al. (1979)] was used. Data for this track is 
given in Appendix A2 . 

For 2D8N modelling, was modified to 5310 MPa and 
average line load intensity is taken as 84.75 kN/m while 
maximum line load intensity is taken as 244.9 kN/m . 

Figures 3.29 to 3.31 show the results of finite element 
analysis using 2D8N routine with the above track parameters for 
Mesh A, Mesh B and Mesh C. Minor principal stress (tr^) shows 
discontinuity at the interface of two layers (Fig. 3.30) 
because o" 3 depends on characteristics of the medium [Burmister 
(1945)]. 

From these figures, it is clear that major (o' 1 ) and minor 
principal stresses have converged for all meshes but 
displacements have converged for mesh B and mesh C only. Hence 
Mesh B may be adopted. 

Mesh A to mesh C were generated by semi automatic mesh 
generation facility. But, even then, it took one to six hours 
to prepare and feed the input data. Hence, it is felt that it 
would be very helpful, if in place of sloping mesh, a 
rectangular mesh which could be generated by automatic mesh 
generation facility can be used. 

In literature, it was found that Selig et al. (1979), and 
Desai and Siriwardane (1982) had used rectangular mesh to model 
the sloping track embankment. Hence, a trial with five 
rectangular mesh types [details are given in Tables 3.4(a) to 
3.4(c) and Fig. 3.32] was carried out with a view to see 




Fig. 3.29 Variation of (7i with depth tor TRACK! 
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Fig. 3.30 Variation of O 3 with depth for TRACK 2 
(2D8N model) 
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Fig. 3.31 Variation of 6 with depth for TRACK 2 
(2D8N model) 
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Table 3.4(a) Meshes used in trials for rectangular mesh selection 


Mesh 

type 

Length of 
Bottom 
Boundary 
(cm) 

Length of 
Lateral 
Boundary 
(cm) 

* 

Total compartments 

Total 
number of 
elements 

Lateral 

Vertical 

D 

455 

400 

7 

8 

56 

E 

665 

585 

9 


90 

F 

665 

585 

11 

14 

154 

G 

1425 

1340 

11 

12 

132 

H 

13700 

13000 

40 

32 

1280 

B 

595 

525 

SLOPING 

120 


* refer Fig. 3.32 









Table 3.4(b) Compartment dimensions of meshes used in trials 
for rectangular mesh selection. 




compartment dimensions in cm 

KSgB 



D 


50.7, 31.7, 31.7, 22.8, 40.6, 76.2, 146 



12.7, 12.7, 12.7, 12.7, 15.2, 35.5, 101, 252 

E 

WLC 

50.7, 31.7, 31.7, 22.8, 40.6, 63.5, 63.5, 127, 153 


WVC 

12.7, 12.7, 12.7, 12.7, 15.2, 22.8, 38.1, 76.0, 152, 



310.0 

F 

WLC 

50.7, 16.5, 15.2, 16.5, 15.2, 22.8, 40.6, 63.5, 63.5, 



127, 153 


WVC 

12.7, 12.7, 12.7, 12.7, 15.2, 12.7, 22.8, 38.1, 38.1, 



63.5, 63.5, 76.0, 101, 183 

G 

WLC 

50.7, 31.7, 31.7, 22.8, 40.6, 63.5, 63.5, 127, 152, 



254, 502 


WVC 

12.7, 12.7, 12.7, 12.7, 15.2, 22.8, 38.1, 76, 152, 



152, 305, 606 

H 

WLC 

FINE TO 


WVC 

COARSE 


* WLC = width of lateral compartment (cm) from left to right 

* wvc = depth of vertical compartment (cm) from top to bottom. 






Table 3.4 (c) Result of trials for rectangular mesh selection 


Mesh type 

'it 

a 3t 

displacement of 


(kPa) 

(kPa) 

top of sleeper 
(cm) 

D 

48.23 

3.07 

3.05 

E 

47.95 

1.82 

3.91 

F 

51.88 

1.65 

4.17 

G 

47.50 

1.46 

5.87 

H 

50.50 

I 

1.37 

10.57 

B 

49.03 

1.77 

3.78 
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Fig. 3.32 Scheme for rectangular mesh discretization. 


Length of the bottom 
boundary 



whether the same results can be obtained for rectangular 
formation as was obtained for sloping formation and to choose 
optimum rectangular mesh size if it gives good results. 

Figure 3.3 3 shows the displacement variation for TRACK 2. 
It seems that as bottom and lateral boundary distances 
increase, displacements also increase. However, mesh E, mesh F 
and mesh B (sloping formation) give same results. Thus , fixing 
bottom and lateral boundary for rectangular mesh around 6.50 
and 6.00 m from center line respectively gives results similar 
to that of sloping formation. 

Figures 3.34 and 3.35 show that results have converged for 
all the meshes including mesh D. Mesh H shows less major 
principal stresses at all depths (Fig. 3.34) as compared to 
other meshes. It is due to the increase of stiffness of ballast 
and subballast layers due to the increase in width of these 
layers for mesh H. 

Thus it can be seen that sloping track geometry can be 
modelled by rectangular mesh without loss of accuracy. Mesh E 
has been chosen as reference mesh for further 2D8N studies 
because, though stresses have converged for mesh D 
displacements are not satisfactory for mesh D. It is shown that 
further refinement does not give better answers and further 
increase in boundaries will give unsatisfactory displacement 
values. The optimum mesh boundaries (mesh E) will give the best 
results. Mesh H (Appendix Bl) which gives most accurate results 
for elastic continuum has been shown to be less accurate for 


track problem. 
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Fig. 3.34 Variation of C7 1 with depth computed by using 
different rectangular mesh discretizationsfor 
TRACK 2 (2D8N model) 
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Minor Principal Stress, 0$ (kPa) 



Fig. 3.35 Variation ot C7 3 with depth computed by 
using ditterent rectangular mesh discre- 
tizations tor TRACK 2 (2D8N model) 



Results obtained by Selig et al. (1979) for FAST TRACK 
section have been reproduced. Measured value of vertical stress 
at top of the subgrade is also shown. 

3.3.5 LIMITATIONS 

2D8N model has following limitations which restrict its 
general purpose use and reduce its versatility : 

1. Effect of rail type on track structure can not be 
evaluated as rail element is absent in the analysis. 

2. In 2D8N modelling, concentrated wheel load has been 
replaced by equivalent line load and the plane strain finite 
element analysis was carried out. Hence the stress distribution 
essentially shows the characteristics of line load. Thus, the 
rapid stress dispersion expected for the concentrated wheel is 
absent for equivalent line load and a very little stress 
dissipation with depth is obtained in this case. This can be 
seen by comparing MULTA results with 2D8N results in Fig. 3.34 
inside the subgrade. It can be seen that 2D8N predicts 
conservative results in terms of vertical stress at deeper 
depths and compares well with measured vertical stresses at the 
top of subgrade in this case. 

This also indicates the necessity of 3 dimensional 
analysis for track problem. However it is obvious that 3-D 
analysis takes more time as compared to 2-D analysis 



[Zienkiewicz and Taylor (1989)] for same problem. Hence an 
attempt has been made to extract 3-D analysis results from 2-D 
analysis . 

At shallower depths, concentrated wheel load gives rise to 
very high stress intensity which cannot be reproduced by 
average line load intensity. Hence, as a modification to 2D8N 
model, following guidelines have been suggested to reasonably 
predict stresses inside the track : 

Major principal stress variation inside the track 
structure due to both average and maximum pressure intensity 
will essentially consist of two parts — upper sharp 
curvilinear part and lower straight vertical part. The 
modification essentially consists of joining roughly the lower 
straight vertical part of average pressure intensity with sharp 
curvilinear part of maximum pressure intensity curve by a 
smooth transition which should be as small as possible (Fig. 
3.36). It is to be noted that average pressure intensity is 
calculated as in Sec. 3.3.3 and can be taken as equal to wheel 
load/spacing between two axles (1.70 m) while the maximum 
pressure is equal to wheel load x 0.43 / width of sleeper. 0.4 3 
indicates the percentage of total wheel load taken by the 
sleeper below the wheel load. 

However, this modification applies to vertical stress 


only. 
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Major Principal Stress, (T, (kPa) 



Fig. 3.36 Proposed modification to 2D8N model. 



3.3.6 COMPARISON OF RESULTS 


Comparison of predicted results using 2D8N model with 
other models and field test results is presented in Sec. 3.5. 

3.3.7 RESULTS AND DISCUSSION 

Parametric studies have been carried out with reference to 
TRACK 1. TRACK 1 is the representation of typical Indian track. 
Data for this track is given in Appendix Al. Results in terms 
of stresses and displacements at important locations are 
presented below for reference track along with the results of 
parametric studies. Analysis has been carried out both with and 
without joint element between different layers. The range of 
parameters used is as follows : 

E fc = 150, 165, 180 MPa 

E sb = 60, 80, 100 MPa 

E = 10, 20, 50 MPa 
s 

d , = 0.15, 0.35, 0.55 m 

sb ' 

Joint element properties used are as given below (also see Sec. 
3.4.6) : 

Joint element property for tie and ballast interface : 

k = 415200 MN/m 3 
n 

k g = 4150 MN/m 3 



Joint element properties for ballast and subballast 
layers : 

k n = 4150 MN/m 3 
k g = 550 MN/m 3 

Joint element properties for subballast and subgrade 
interface : 

k n = 550 MN/m 3 
k s = 165 MN/m 3 

E^. 1 for discrete real system was changed to E t2 for an 
equivalent continuous system as given below (Sec. 3.3.3) 

E t2 = °’ 249 x 10550 / °- 65 
= 4041.5 MPa 


Wheel load was converted into equivalent line load as follows : 
Average line load = 159.4/1.7 

= 93.77 KN/m 

Maximum line load = 0.43 x 159.4 / 0.249 

= 275.27 KN/m 

Normal displacement was restricted on sides and bottom while 
top edge was kept free. Two bottom corners were assumed as 
fixed. 

Mesh E was used for discretization, in general. 

Results of parametric studies have been tabulated in Table 3.5. 
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Figure 3.37 shows displacement under the load with depth 
for the reference track by 2D8N with and without joint 
elements. Displacement remains fairly constant till top of 
subgrade as the tie, ballast and subballast layers have very 
low compressibility. 

It has been noted that vertical and horizontal stress 
values at the top of subgrade remain same as that of major and 
minor principal stress values respectively. Hence, for 2D8N and 
3D20N model, only variation of major and minor principal stress 
values has been presented. 

Figures 3.38(a) to 3.38(c) show distribution of major and 
minor principal stresses and maximum shear stress with depth 
for reference track for 2D8N model with and without joint 
elements. Major principal stress predicted by modified 2D8N 
model is also shown. Major principal stress distribution is a 
smooth curve through all the layers. However, minor principal 
stress depends on media characteristics while major principal 
stress does not [Burmister (1945)]. Within one layer, minor 
principal stress decreases from top to bottom. Minor principal 
stress is tensile in ballast and subballast layers. While, in 
subballast layer, small tension (up to 50 KPa) may be allowed 
but ballast layer being fully granular cannot sustain any 
tension. This is the disadvantage of elastic analysis. This 
also calls for elasto-plastic analysis. However, later in Sec. 
3.5, it is shown that results predicted from elastic analysis 
compare well with results from elasto-plastic analysis inside 
the subgrade. As the maximum shear stress depends on minor and 
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Displacement , S (mm) 



Fig. 3.37 Variation of displacement with depth for 
TRACK 1 (2D8N model) 
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Major Principal Stress, Qi (kPa) 



Fig. 3.38(a) Variation of CJi with depth for TRACK 1 
(2D8N and modified 2D8N model ) 
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Minor Principal Stress, Cr 3 (kPa) 




Fig. 3.38(b) Variation of (J 3 with depth for TRACK 1 
(2D8N model) 
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^max (kPa) 



Fig, 3.38(c) Variation of *C m ax with depth for TRACK 1 
(2D8N model) 






major principal stresses, maximum shear stress curve also has 
discontinuity at the interfaces. 

From Figures 3.38(a), it may be seen that the dispersion 
of major principal stress with the depth is poor. Stresses 
decrease by only 20 percent for 1 m depth inside the subgrade. 
This is because of the assumption of line load in this model. 

Figure 3.39 shows major and minor principal stresses and 
Shear stress distribution on top of the subgrade with the 
length in x-direction (but inside the subgrade) . It may be seen 
that maximum major principal stress occurs under the rail seat. 
As major principal stress is equal to vertical pressure on top 
of subgrade, maximum vertical pressure on top of subgrade is 
also located under the rail seat. 

Figure 3.40 shows major and minor principal stress on top 
of subgrade under the rail seat (cr^ and cr^ respectively) for 
different soil types, for three different values of thickness 
of subballast layer. Thickness of ballast layer is kept 
constant (0.25 m) in all three cases. It can be seen that cr t 
increases with increase in Young's modulus of subgrade, cr^ 
decreases with increase in Young's modulus of subgrade. 
However, both stresses decrease with increase in subballast 
depth. Results have been plotted with and without joint 
elements. The effect of joint elements is to reduce both types 
of stresses. 

Figure 3.41 shows the displacement of top of sleeper for 
different soil types and for different depths of construction 
with and without joint elements. It can be seen that depth of 



Distance (in x-direction) (m) 
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Fig. 3.39 Variation of major and minor principal stress along the sleeper 
on top of subgrade for TRACK 1 (2D8N model) 
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Young's Modulus of Subgrade, E s (MPa) 

Fig. 3.40 Variation of (In and Q 3t with E s for diff- 
erent depths of subballast (2D8N model) 
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Young’s Modulus of Subgrade, Eg (MPa) 

Fig.3.41 Variation of displacement on top of sleeper 
with E s for different depths of subballast 
(2D8N model) 
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construction has very little effect on displacement whereas 
type of soil has very significant effect on displacement of top 
of sleeper. Effect of joint elements is to increase the 
displacement. 

Figure 3.42 shows that the gain achieved by increasing the 
depth of ballast layer but keeping total depth of formation 
constant, in terms of reduction of settlement or major 
principal stress on top of subgrade is not significant. 



C r n (k Pa) 
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8N model) . 


Displacement on top of sleeper 
8 (mm) 



3.4 THREE DIMENSIONAL FINITE ELEMENT ANALYSIS 


3.4.1 INTRODUCTION 

Three dimensional finite element analysis of track has 
been carried out in order to calculate deflections and stresses 
inside the subgrade. Railway track model for three dimensional 
finite element analysis has been described below in detail. 
Stress distribution and displacements for a typical track have 
been presented and the parametric studies are carried out using 
three dimensional finite element analysis with reference to 
this track. 


3.4.2 DESCRIPTION OF THE THREE DIMENSIONAL FINITE ELEMENT 
ROUTINES 

Following 3-D finite element routines have been developed 


1. 3D8N 

2. 3D20N 

Both routines are based on displacement formulation of 
elasticity eguations. Both routines use linear, elastic 
material properties. 3D8N routine uses 8 noded linear 
isoparametric brick elements. 3D20N routine uses 20 noded 
isoparametric brick elements. 

Both routines were tested successfully for standard 



Boussinesq's problem and standard cantilever beam problem. 
Comparison of results from 3D20N routine (which has been used 
in this analysis) for Boussinesq's problem is given in Appendix 
B2 and comparison with exact solution of cantilever beam 
problem is given below. Comparison of results with exact 
solution has shown that 3D20N gives good results with quite 
less number of elements while 3D8N gives unsatisfactory results 
even with very high number of elements. 

For example, consider the following standard cantilever 
beam problem [Fenner (1986)] : 

Input data : 

Cantilever beam 

Deflection calculated at outer tip of beam 

■v = °- 3 

E. = 1000000 
bm 

Table 3.6 shows the comparison of results from 3D8N vis - a - 
vis 3D20N model for this problem. 

It is clear from Table 3.6 that even with 32 elements, 
3D8N is not able to give good results because it has linear 
shape functions which cannot model bending satisfactorily while 
3D20N routine gives nearly exact answer with only 2 elements. 
Also, same observations were noted for standard Boussinesq's 
problem [Yudhbir et al . (1991)]. Hence, 3D20N routine has been 

used for further finite element analysis. 

3D20N routine can handle layering of different materials 



Table 3.6 Beam problem 3D8N vis-a-vis 3D20N 


Program used 

Mesh refinement 

FEM answer 

Exact answer 

3D8N 

2 elements 

0.00190 

0.00864 


4 elements 

0. 00435 



“w °- 15 > 

(0.00439) 



8 elements 

0.00651 



16 elements 

0. 00745 



32 elements 

0.00772 


3D2 0N 

2 elements 

0.00860 



* Above calculations are independent of unit system 













and all types of boundary conditions. The routine has facility 
to simulate joint between different layers by use of joint 
element of zero thickness [Goodman et al. (1968)]. 16 noded 
surface joint element based on relative displacement 
formulation has been used [Burgohain and Shah (1978)]. Joint 
element subroutine was checked for the algorithm given by 
Burgohain and Shah (1978) . 

3D20N routine has capability to analyze the soil structure 
interaction problems. It has facility to combine beam element 
with soil element. This routine has been checked by comparing 
results with MULTA model [Selig et al. (1979)] which also 
considers beam as 1-D element (Figs. 3.61 and 3.62). 

Stress extrapolation from Gauss points to nodal points has 
been carried out by bilinear extrapolation. 

Three types of data input facilities are available. In 
fully manual data type, data related to nodal point 
coordinates, element connectivity, etc. has to be fed for each 
and every point. But, main advantage of this routine is that 
one can prepare data for virtually any type of mesh 
discretization. In semiautomatic data input, data related to 
some nodal points and elements is needed. This can generate 
only linear geometric mesh patterns. Fully automatic data input 
facility can only be used to generate rectangular meshes. Data 
input consists of basic parameters (e.g. total number of nodal 
points, material types, elements etc.), material properties, 
total number of compartments in x, y and z direction and length 
of each compartment, boundary conditions on six faces and 8 




3. 4. 3. 2 Boundary conditions 


The boundary conditions adopted for 3-D analysis were 
similar to that of 2-D analysis and are given below : 

At all faces, displacements in the normal direction were 
restricted. 

At all corner lines, displacements across the line were 
restricted. However, displacements in the direction of line 
were allowed. 

At all corner nodes (except for top center node) , all 
displacements were restricted. 

On top center node, only vertical displacement was 
allowed. 

3. 4. 3. 3 Rail element 

Rail has been replaced by one dimensional beam element 
having two degrees of freedom (displacement in vertical 
direction and rotation) at each node for finite element 
analysis. 1-D beam element has advantage over 3-D rail section 
of equivalent El [Desai and Siriwardane (1982)] because if rail 
is modelled by 3-D element of equivalent El then assumption 
will have to be made regarding contact area between rail and 
sleeper. As rail is supported directly over the sleeper for 
considerable distance, this area presumably, may have important 
bearing on stresses transferred. Thus, ratio of height to width 
of rectangular rail section chosen will be empirical and may 
affect direct stresses transferred to the sleepers. It can lead 
to underprediction of stresses and displacements at the top of 



subgrade (Figs. 3.44 and 3.45). 

In Figs. 3.44 and 3.45, displacement and stress 
distribution under the rail seat along the depth are shown. Two 
cases as given below are considered and in each case, El of 
rail section is made equal to actual El : 

1. 1-D rail section 

2. 3-D square rail section. 

Also shown are the results predicted by MULTA and 3- 
dimensional finite element analysis by Desai and Siriwardane 
(1982) . Desai and Siriwardane (1982) results and 1-D rail 
results are different due to the assumptions about rail 
section. However, 1-D rail gives best results compared to MULTA 
model which is like a closed form solution. 

1-D rail also reduces computational time by reducing the 
total number of nodal points and the band width of global 
stiffness matrix. 

Tie and rail connection for 1-D rail has been arrived at 
as follows : 

Two cases have been considered regarding assumption of tie 
and 1-D rail connection : 

a. Node where tie is connected with rail is not restrained 
in any direction. 

b. Node where tie is connected with rail is restrained in 


lateral direction. 
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Fig. 3.44 Displacement under the rail seat versus 
depth for different rail conditions for 


TRACK 2 
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Fig. 3.45 0“i versus depth for different rail conditions 
for TRACK 2 



It has been assumed that rail is free to rotate and 


displace vertically at the connection between rail and tie in 
both cases. Results for both these cases are shown in Figs. 
3.44 and 3.45. Results also show stress distribution without 
rail element. Thus, the effect of presence of rail element can 
be seen. For further analysis, it is assumed that the rail is 
restrained in lateral direction. 

3. 4. 3. 4 Sleepers 

Sleepers have been modelled by 3-D rectangular elements. 
Discontinuous nature of sleeper could be maintained here. 
However, to maintain the rectangular geometry, interspaces 
between adjacent sleepers (cess) and spaces on side of sleepers 
(crib) were filled with material of very low compressibility 
(E f = 0.7 MPa) . 

A detailed parametric study was carried out to calculate 
E^. The effect of filling cess and crib will obviously increase 
the stiffness of track and hence stresses on top of subgrade 
will reduce. Hence, a parametric study (using TRACK 2, single 
wheel load ; no joint element) was carried out with E^.= 0.007, 
0.7, 210 MPa. The corresponding values of cr^ were obtained as 
44.57, 44.92 and 41.62 kPa. E f = 0.7 MPa gave slightly higher 
results as compared to E f = 0.007 MPa due to rounding off 
errors during calculations. Thus, E^ = 0.7 MPa was adopted as 
the Young's modulus of filling material. 
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3.4.4 SELECTION OF FINITE ELEMENT MESH 

Rectangular track outline is adopted for 3-D analysis. The 
symmetry in x and y directions is made use of and hence only 
one quarter of track section need to be considered. 

For choosing optimum mesh size, five types of mesh 
discretizations were studied. Tables 3.7(a) to 3.7(c) show the 
properties of five meshes used and the results obtained. Table 
3.7(c) also shows the CPU time taken by these meshes to analyze 
TRACK 2 with single wheel load. Joint elements are not used. 
Mesh B' is recommended for finite element analysis. However, 
this mesh does not give sufficient points after 1 m depth to 
give clear picture of stress distribution. For this, mesh C' is 
recommended. Mesh C' is shown in Fig. 3.46. 

It can be seen that results approach towards MULTA results 
which is a type of closed form solution but never underpredicts 
the stresses as in the case of Desai 's 3DFEM model [Desai and 
Siriwardane (1982)]. In assuming 1-D rail , stresses dissipated 
by rail element by direct shear mechanism have been neglected 
and only bending resistance of rail element is considered. This 
gives safer and better results as compared to substituting rail 
by rectangular element of equivalent El [Desai and Siriwardane 
(1982) ] with some assumption about contact area between sleeper 
and rail and thickness of rail element. 



Table 3.7(a) Meshes used in trials for 3-D mesh selection 


Mesh type 

x x y x z 

x-boundary 

(cm) 

y-boundary 

(cm) 

z -boundary 
(cm) 

A' 

3 x 4 x 5 

442 

320 

381 

B' 

4 x 6 x 6 

442 

320 

381 

i 

C' 

5 x 6 x 8 

442 

320 

381 

D' 

7 x 8 x 10 

660 

420 

650 

E' 

4 x 6 x 6 

240 

190 

275 


*xxyxz= x, y, and z denote total number of compartments 

in x- , y- and z- direction respectively. 







Table 3.7(b) Compartment dimensions for meshes used in trials for 
3-D mesh selection. 


Mesh Compartment Compartment dimensions in cm 
type name 

X-COMPARTMENT 82.5, 54.5, 305.0 

A' Y-COMPARTMENT 12.7, 25.4, 25.4, 254 

Z-COMPARTMENT 17.8, 38.1, 15.2, 50.8, 254 

X-COMPARTMENT 41.9, 40.6, 54.5, 305,0 

B' Y-COMPARTMENT 12.7, 25.4, 25.4, 25.4, 25.4, 203.0 

Z-COMPARTMENT 17.8, 17.8, 20.3, 15.2, 50.8, 254.0 

X-COMPARTMENT 66.0, 16.5, 16.5, 38.0, 304.8 

C' Y-COMPARTMENT 12.7, 25.4, 25.4, 25.4, 25.4, 203.0 

Z-COMPARTMENT 17.8, 17.8, 20.3, 15.2, 15.2, 50.8, 101, 152 

X-COMPARTMENT 33.0, 33.0, 16.5, 16.5, 38.0, 114.0, 406.0 
D' Y-COMPARTMENT 12.7, 25.4, 25.4, 25.4, 25.4, 25.4, 25.4, 254 

Z-COMPARTMENT 5.08, 12.7, 17.8, 20.3, 15.3, 15.3, 30.5, 76, 

152, 306 

X-COMPARTMENT 41.9, 40.6, 54.5, 101 

E' Y-COMPARTMENT 12.7, 25.4, 25.4, 25.4, 25.4, 76.2 

Z-COMPARTMENT 17.8, 17.8, 20.3, 15.2, 50.8, 152 













Table 3.7(c) Result of trials for 3-D mesh selection 


Mesh type 

CPU time taken 
on convex 220 

ff lt 

Rail displacement 


(sec) 

(kPa) 

(cm) 

A' 

74.6 

56.0 

0.240 

B' 

437.0 

46.5 

0.225 

C' 

1774.5 

46.0 

0.220 

D' 

12192.0 

45.0 

0.250 

E' 

433.0 

48.0 

0.195 










1D-Rail element 



1 cm = 50 cm 


Fig. 3.46 Mesh C* discretization. 




3.4.5 LIMITATIONS 


Limitations of 3D20N model are discussed below : 

1. 3D20N model is a linear, elastic 3-D finite element 
model and predicts stresses and displacements inside the 
subgrade. 

Selig et al. (1979) have carried out the dynamic testing 
of track and found that one of the significant features shown 
by dynamic records is recoverable deformations of the track 
system under instantaneous loads and recommended the use of 
elastic models for predicting the dynamic response of the track 
system. They have also noted that the response of the track 
system is approximately linear. Hence assumption of linear 
elastic model may not be the limitation. 

2. However, Janin et al. (1983b) pointed out that for 
railway track where a material of high modulus overlies one 
with a lower modulus, tensile horizontal stresses develop in 
the bottom of the top layers. One of the limitations of the 
elastic analysis is that the materials considered can support 
tension as well as compression. Thus, elastic analysis predicts 
tensile horizontal stresses in the granular layer which is 
clearly unrealistic. Hence, simple elastic analysis is 
unsuitable for predicting horizontal stresses in top layers. 

However, Janin et al. (1983b) also conclude that elastic 
analysis is satisfactory for predicting vertical stresses 
acting on the subgrade. In Sec 3.5, it is shown that the 



horizontal stresses predicted by elastic and elasto-plastic 
analysis inside the subgrade are in the same range. Moreover, 
it is to be noted that elasto-plastic analysis takes 10 times 
more time than elastic analysis. 

Hence, unless the real mechanism of subgrade failure is 
understood (after close monitoring of track operation on 
subgrade for few years) , elastic analysis may be the better and 
practical choice. 


3.4.6 EVALUATION OF PARAMETERS USED IN CONSTITUTIVE 
RELATIONSHIP 

3D20N analysis uses linear elastic material properties and 
hence Youiig's modulus and Poisson's ratio are the parameters to 
be evaluated for different materials. 

In this study, two tracks have been considered. For TRACK 
2, parameters suggested by Selig et al. (1979) have been used. 

TRACK 1 is a typical Indian track. The best procedure for 
input parameters would have been to study constitutive 
modelling for media in track support structures as done by 
Janardhanam [Desai and Siriwardane (1982)]. However, in the 
absence of such data, ranges of parameters suggested by RDSO 
through personal discussion have been used. 

To determine the Young's modulus of subgrade, cyclic plate 
load tests and cross hole wave propagation tests were carried 
out by RDSO. The detailed results are reported in Kumar et al. 
(1990), Kumar et al. (1991a) and Kumar et al. (1991b). Figure 



3.47 shows the results of cross hole wave propagation tests on 
the running track for Sandilla-Rahimabad section near Lucknow 
at five different locations. However, it was noted that at 
short depths below formation level, soil might be mixed with 
ballast and moorum and hence, at these depths, Young's modulus 
of subgrade measured is higher. Cyclic plate load test was 
carried out on top of formation and hence reflects the same 
trend (Table 3.8). Hence, Young's modulus of subgrade was taken 
as 20 MPa for the reference track. 

Joint elements have been used between tie-ballast, 
ballast-subballast and subballast-subgrade interfaces. The 
values of joint stiffness have been determined as follows : 

The following guidelines can be noted in literature. 

1. A high joint stiffness is a value which yields 
negligible joint displacements when compared to the elastic 
displacements of the layers; a low stiffness is a value which 
leads to joint displacements that dwarf the elastic 
displacements of the blocks; a moderate joint stiffness 
corresponds to displacements of joints of the same order of 
magnitude as the elastic displacements of the layers [Goodman 
et al. (1968) ] . 

5 

2. Value of k n is normally 10 times the maximum value of 
Young's modulus of the bodies in contact to avoid significant 
overlapping. For 'no slip' condition, k gl and k s2 are given 
high values similar to k n and for perfectly smooth contact with 



Young's Modulus of Subgrade, E$ (MPa) 



Fig. 3.47 Es versus depth below formation level -wave propagation 
test results on Sandila - Rahimabad section. 
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Table 3.8 Summary of cyclic plate load test results (Tests 

carried out on Sandila-Rahimabad section near Lucknow) 



1114/3 53.5 
1114/16+2 38.5 
1114/17 55.2 
1116/8+9 62.0 


* Location indicates section on the operating track. 





no resistance to slip, k = k „ = 0. In case of separation, k 

s ± s /c n 

= k gl = k s2 = 0 [Burgohain and shah (1978)]. 


3. It has been observed from the examples given in 

literature that moderate values of k are of the order of E/1 

where E is Young's modulus of adjoining layer and 1 is the 

thickness of the block while low values of k , and k _ are of 

si s2 

the order of 1/10 th of k n (moderate) [Goodman et al . (1968)]. 


From personal discussion with RDSO, the following range of 
joint stiffnesses has been recommended : 


Interface 

tie-ballast interface (1) 

ballast-subballast 
interface (2) 

subballast-subgrade 
interface (3) 


high 


moderate 


moderate 


si 

small 

moderate 


moderate 


small 


moderate 


moderate 


Next, a detailed parametric study was carried out to 
evaluate the effect of each joint stiffness. Following 
procedure was adopted for parametric study : 

With reference to TRACK 2, according to guidelines given 
above the following track parameters were taken as reference 


values : 



Interface 

k n (MN/m 3 ) 

k sl (MN/m 3 ) 

1 

415000 

4150 

2 

8300 

900 

3 

1400 

550 

k , 

= k „ 


si 

s2 


The range 

of values varied 

is as follows 


Interface 1 

k R (MN/m 3 ) 41500000., 415000., 8300., no joint element 

k gl (MN/m 3 ) 4150., 550. 

Interface 2 

k n (MN/m 3 ) 8300., 550. 
k gl (MN/m 3 ) 900., 280., 55. 

Interface 3 

k n (MN/m 3 ) 1400., 550., 140. 

k gl (MN/m 3 ) 1400., 550., 55. 

The above ranges are chosen because for very high values, 
results become constant and for very low values instability 
occurs (i.e. change in displacement becomes significant) . 

Keeping in mind that displacement of tie, ballast and 
subballast layers is of the order of 0.2 to 0.4 mm and that of 
subgrade is 1.55 mm and following the guidelines given by RDSO 



and using results of parametric study, the following values of 
joint stiffnesses were arrived at : 




k n (MN/m 3 ) 

k sl (MN/m 3 ) 

k s2 (MN/m 3 ) 

Interface 

1 

415000. 

4150. 

4150. 

Interface 

2 

4150. 

550. 

550. 

Interface 

3 

550. 

165. 

165. 


3.4.7 COMPARISON OF RESULTS 

Comparison of predicted results using 3D20N model with 
other models and field test results is presented in Sec. 3.5. 


3.4.8 PARAMETRIC STUDIES 

Parametric studies have been carried out for a range of 
parameters with reference to TRACK 1. TRACK 1 is a 
representation of the typical Indian track. Details of this 
track given in Appendix A1 have been arrived at after personal 
discussion with RDSO and from published RDSO reports. 

Parametric study has been carried out with reference to 
this track by changing the following track parameters. 

Young's modulus of ballast = 150, 180 MPa 

Young's modulus of subballast = 60, 80, 100 MPa 

Depth of subballast = 0.15, 0.35, 0.55, 0.75 m 


Young's modulus of soil =10, 20, 30, 50 MPa 



. 4 

Moment of inertia of rail = 1600, 2158, 3055 cm 

Tie spacing = 0.50, 0.65, 0.76 m 

Making use of symmetry, only one-fourth of the track 
section has been discretized. Mesh C' shown in Fig. 3.4 6 has 
been used to model the reference track. For various sub-ballast 
thicknesses, mesh division in z-direction was varied. Rail 
element has been modelled as one dimensional beam with 
equivalent El as that of rail section. Joint elements were used 
at all the three interfaces viz. tie-ballast, 
ballast-subballast and subballast-subgrade. The properties of 
joint elements used are : 



k n 

k sl 

k s2 

(all in 

Tie ballast interface 

415000. 

4150. 

415. 

MN/m 3 ) 

Ballast subballast interface 

4150. 

550. 

55. 


Subballast subgrade interface 

550. 

165. 

165. 



Boundary conditions used are same as mentioned in Sec. 
3. 4. 3. 2. 

Double wheel load has been considered throughout this 
analysis . 

3.4.9 RESULTS AND DISCUSSION 

In this section, first the trends observed for reference 
track are presented and then the results of parametric studies 



have been discussed. 


Figures 3.48 to 3.52 show the typical results in terms of 
stresses and displacements for reference track, TRACK 1. 
Figures 3.48 to 3.51 show results for three cases viz. double 
wheel load with joint elements, double wheel load without joint 
elements and single wheel load without joint elements. This 
brings out the effect of double wheel load and joint elements. 

It was found that due to the vertical nature of wheel 
load, vertical stresses and major principal stresses as well as 
lateral stresses and minor principal stresses are same for 
single wheel load case. Hence, in this study, only principal 
stresses are considered. 

Figure 3.48 shows the displacement under the load. Figures 
3.49(a) and 3.49(b) show the major and minor principal stresses 
inside the subgrade. Minor principal stress distribution has 
been shown inside the formation and subgrade. As granular 
material cannot sustain tension, minor principal stress in 
granular material has no meaning. It has been shown in Sec. 3.5 
that stresses inside the subgrade calculated by elastic 
analysis match well with the stresses calculated by 
elasto-plastic analysis which takes into consideration 'no 
tension ' condition in granular material. 

It may be seen (Fig. 3.49a) that for single wheel load 
case, cr 1 decreases from 53 kPa on top of subgrade to 17 kPa at 
1.0 m inside the subgrade. This shows that 68 % stress 
dispersion takes place within 1 m inside the subgrade while for 
double wheel load, it is only 14 % (from 42 kPa to 36 kPa) . 
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Displacement, 6 (mm) 

0 2 4 6 8 10 12 



Fig. 3.48 Variation of displacement under the rail seat 
with depth tor TRACK 1 (3D20N model) 
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Major Principal Stress, dj (kPa) 


0 40 80 120 160 200 240 280 



Fig. 3.49^Variation of (7i with depth for TRACK 1 
(3D20N model) 



Minor Principal Stress, CT 3 (kPa) 




Double wheel load with joint element 
Double wheel load without joint element 
Single wheel load without joint element 


Fig. 3.49(b) Variation of O 3 with depth for TRACK 1 
(3D20N model) 



165 


16.25 1 


A 


zz: 


2 


Depth from top ol 

(< 

l 

sleeper 

:m) 

\ i 

-JL _ 8 

f 

KOTSG)_ 


o 

Q- 

in 

i/> 

a> 

«*■» 
in 

o 

a 

• 

o 

.E 160 


N X 


120 


o 

o' 


200 


240 


280 


\ <x \ 46(OTSB >| 


8K0T5G) 

\ \v \ 

80 -\\\ 

A \ 

\\ ~ ... 

\ \ s 46 <°^l 

\ 

\ 


\ 




^y I 

v 25 (OTB) / 

\ / 

\ / 

X / 

\ / 

V ' 

2 5 (OTB) 

—• — Double wheel load without joint element 
— — — Single n n n // n 


OTSG- On top of Subgrode 
OTSB - On top of Subballast 
OTB - On top of Ballast 


_L 


± 


0 40 80 120 

Length Along Sleeper 


Fig. 3.50 Variation of major principal stress 
along the sleeper at different depths 
for TRACK 1(3D20N model) 
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Fig. 3.51 Variation ot rail deflection along the rail direction for TRACK 
(3D20N model) 
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Fig. 3.52 Stresses on top of subgrade along the 
sleeper for single wheel load with joint 
element (3D20N model ) 




It may be seen (Figs. 3.48 and 3.49) that effect of joint 
element is to reduce stresses and to increase displacements in 
top layers. Joint element reduces by 8 % on top of subgrade. 
Joint element reduces a 3 from 10 kPa to 7 kPa on top of 
subgrade while 5 on top of sleeper increases from 5.8 mm to 6.2 
mm (increase by 7 %) . Desai and Siriwardane (1982) have also 
observed increase in 6 on top of sleeper of similar magnitude 
(for 147.15 kN load, 6 increases by 6 %) . 

Effect of double wheel load is to reduce the stresses at 
shallow depths and increase the stresses at deeper depths 
(below 0.80 m in this case). Double wheel load increases 6 on 
top of sleeper from 3.9 to 5.8 mm (i.e. by 8 %) while cr^ 

increases from 5 to 10 kPa on top of subgrade, cr^ decreases 
from 248 kPa to 178 kPa on top of ballast for double wheel load 
case but the changes are negligible on top of subgrade. 

Figure 3.50 shows major principal stress distribution 
along the sleeper at depth egual to 0.46 m and 0.81 m. Maximum 
stress on top of subballast and ballast is found to be slightly 
on the right side of the rail seat. However, on top of 
subgrade, maximum major principal stress occurs below the rail 
seat. These findings are similar to those observed by Selig et 
al. (1979). 

Figure 3.51 shows the rail deflection along the rail (in 

the y-direction) . Figure 3.52 shows the distribution of major, 

intermediate (<x_) , minor principal stresses and maximum shear 

stress (t ) on top of subgrade for single wheel load case 
max 

with joint elements. It may be noted that maximum major 



1 ' 


principal stress occurs below the rail seat but r _ on top of 

max 

subgrade is maximum slightly to the right of the rail seat. 

However, the difference between maximum x and t _ below the 

max max 

rail seat is insignificant and for all practical purposes, T max 
below the rail seat may be taken as maximum T max on top of 
subgrade . 

As already mentioned above, four types of track structure 
and four types of subgrade were selected for parametric 
studies. Each track structure was put on each type of subgrade 
and 3-D finite element analysis was carried out. Track 
structure 2 with subgrade 2 (TRACK 1) has been used as 
reference track for parametric studies. These studies were 
carried out to see the effect of tie spacing, moment of inertia 
of rail, quality of ballast and quality of subballast on the 
major principal stress on top of subgrade. Results of 
parametric studies are given in Table 3.9. 

A non-dimensional parameter 'm ' has been introduced as a 
measure of ballast , subballast and subgrade characteristics as 

E 

equ 

m 

E s 

where 

E = equivalent Young / s modulus of ballast 
subballast layer and is expressed as 


(3.20) 


and 
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Table 3 9 CContd D Results of parametric studies C3D20N model D 














(3.21) 


E b d b + E sb d sb 


d.+d , 
b sb 


with 


= Young's modulus of ballast, 

= Young's modulus of subballast, 

= thickness of ballast layer, and 
= thickness of subballast layer. 


E can be derived on basis of stress distribution in the 

equ 

track system as follows : 

Figures 3.53 and 3.54 show the assumed stress distribution 
in the track system for actual and equivalent system. 
Therefore, 


( d b + d sb ) tan 0 = d b tan e i + d sb tan e : 


(3.22) 


Let us assume. 


d b E b + d sb E sb = < d b + d sb> E 


(3.23) 


E b d b + E sb d sb 
d b + d sb 


(3.24) 


From Eqs. (3.22) and (3.24), 


E b d b + E sb d sb 
d b tan 8 1 + d sb tan ®2 


E / tan 6 - 



SLEEPER 



Fig, 3.53 STRESS DISTRIBUTION IN ACTUAL TRACK SYSTEM 
(TADATOSHI ET AL V A> 



Fig. 3.54 Assumed stress distribution in equivalent 


track system. 



(d fa tan © 1 / d sb tan © 2 + l) tan © 2 


(3.25) 


For a given track structure, 0^ and © 2 are constant. 

Hence E oc tan 0 (3.26) 

^egu can also 3 be obtained by assuming track consisting of 
composite elastic plates as shown below : 


For composite plates, neutral surface is located by (Fig. 


3.55) 


y n = 


E 1 ?! A 1 + E 2 y 2 A 2 
E 1 A 1 + E 2 A 2 


(3.27) 


•• (E-l + e 2 A 2 )y n = E 1 Y 1 A 1 + e 2 y 2 a 2 


(3.28) 


Thus E A = E. A. + E_ A_ 

egu equ 11 2 2 


(3.29) 


If b is constant, 


and 


d egu - d l + d 2 ' 


E 1 d 1 b + E 2 d 2 b = E d b 


" E 1 d l + E 2 d 2 * E 4 



Plate 

1 

— Ballast 

Plate 

2 

— Subballast 



Fig, 3.55 ELASTIC PLATES 



_ E 1 *1 + E 2 *2 


(3.30) 


equ 


d l + d 2 


This is same as Eq. (3.21). However this assumption is 
good for elastic analysis only because it enables the plate to 
take tension. 

Results of parametric studies are shown in Figs. 3.56 to 
3.60. Figures 3.56 and 3.57 show the effect of Young's modulus 
of subballast layer on top subgrade stresses. Figure 3.57 shows 
the effect of Young's modulus of ballast layer on stresses at 
the top of subgrade. Increase in E gfc from 60 MPa to 100 MPa 
decreases tr lt by 2 % while increase in E^ from 150 MPa to 180 
MPa decreases <r ^ by only 1 %. 

Figure 3 . 58 shows variation of top subgrade stresses with 

moment of inertia of rail. Increase in moment of inertia of 

4 ... 

rail from 1600 to 3055 cm decreases cr lb by 5 %. This finding 

is same as reported by So (1978). Figure 3.59 shows the 

variation of top subgrade stresses with tie spacing. Tie 

spacing has moderate effect on stresses on top of subgrade 

(about 6 to 17 percent for change in tie spacing from 0.50 to 

0.76 m) . 

Figure 3.60 shows the variation of stresses on top of 
subgrade with m for four depths of construction. As m 
increases, major principal stress on top of subgrade decreases. 
Also, as depth of construction increases, major principal 
stress on top of subgrade decreases as expected. As noted in 
Sec. 2.4, there are two aspects in track stability namely the 



Olt < kp a) <5' cr 1t (kPa) 
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E s b (MPa) 

3.56 Variation of On with Young's modulus of 
subballast, £ sb for reference Track (3D20N) 



E b (MPa) 

Fig. 3.57 Variation of Oit with Youngs modulus of 
Ballast, E b for reference Track (3D20N) 






Fig. 3.59 Variation of Oit with tie spacing for TRACK 1 
(3D20N) 
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Ffg. 3.60 Variation of ( Tit with 



track equipment which consists of rail and tie spacing and the 
formation quality which is related to quality and thickness of 
better materials and the quality of subgrade. It can be 
observed that if m is changed by changing the formation quality 
either by changing quality or thickness of formation or quality 
of subgrade, it will follow the trend shown in Fig. 3.60 for a 
particular track equipment. However, if track equipment is 
changed, stresses will not follow the trend shown there. Thus, 
in Fig. 3.60, change in major principal stress on top of 
subgrade due to the change in Young's modulus of subballast 
layers will follow the curve for that particular depth. 
However, if tie spacing is changed, then major principal stress 
observed, will not follow the trend in Fig. 3.60. Thus 
formation quality is represented by m. 

Figure 3.60 also shows the results from Janin et al. 
(1983c) predicted by elasto-plastic analysis by Rosalie model. 
These results also show straight line relations. Two sets of 
results can be compared by making following corrections. 

a. Magnitude of load : Janin et al. (1983c) presented 
results for wheel load = 10 t. 

b. Janin et al. (1983c) carried out calculations for 
single wheel load. 

However, discrepancy in results occurs which may be due to the 
following reasons : 



1. For results given by Janin et al. (1983c) , subgrade 
depth has been taken as only 0.15 m thick and bottom is assumed 
to be restrained. 

2. Janin et al. (1983c) presented results for track 
consisting of four different layers with an additional gravel 
layer. However, in the absence of sufficient data about 
thickness of different layers, an average value of E e ^ u = 150 
MPa has been assumed. In reality, E e ^ u will be changing with 
the structure type. 



3.5 COMPARISON OF RESULTS 


The three models described above viz. ELASTIC. F, 2D8N 
(also modified form) and 3D20N have been compared with 
published results predicted by other models and measured in the 
field. Results of four track studies viz. FAST TRACK (TRACK 2), 
So (1978) track, Vienna Arsenal single sleeper and Derby panel 
full scale model tests and field measurements carried out by 
Research, design and specifications organization, Indian 
railways (RDSO) have been compared. As far as possible, 
comparison of results has been made for displacements, major 
and minor principal stresses. For TRACK 1, results predicted by 
using the three models have been compared with results obtained 
by published empirical formulae. For comparison, 4x6x6, 4 x 
6x8 and 4 x 6 x 10 mesh discretizations have been used. 

3.5.1 FAST TRACK 

Selig et al. (1979), and Desai and Siriwardane (1982) have 
reported stresses and displacements using elastic model based 
on Burmister's multilayer theory and finite element analysis 
respectively on FAST TRACK. The track description for FAST 
TRACK is same as for TRACK 2 and is given in Appendix A2 . 

Figures 3.61 and 3.62 show vertical displacement and 
vertical stress distribution below the wheel load along the 
depth for FAST TRACK using various models mentioned therein. 
Using the same track parameters, results were predicted for 
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single wheel load using ELASTIC. F, 2D8N and 3D20N models. Joint 
elements were not used. It is clear that 3D20N predicts results 
similar to MULTA model which is like a closed form solution. 
Stress prediction by using ELASTIC. F model is also very close 
to MULTA model. However, it can be seen that 2D8N model 
predicts lower stresses at shallower depths. Modified 2D8N 
gives results that are quite satisfactory. 

Figure 3.63 shows the stress distribution along the 
sleeper by using the various models. At subgrade top (71 cm 
depth), both 2D8N and 3D20N give results close to MULTA model. 
However, at shallower depths (0.18 m) , 2D8N predicts lower 
stresses. 3D20N gives results similar to MULTA for all depths. 
Figure 3.64 compares displacement along the rail predicted by 
3D20N model with MULTA, PSA and ILLITRACK models. 

3.5.2 SO (1978) TRACK 

So (1978) has carried out finite element studies and has 
predicted results in terms of variation of maximum vertical 
stress on subgrade for change in tie spacing, ballast depth and 
young's modulus of subgrade with reference to the typical 
track. 

The properties of reference track along with the ranges of 
parameters used by So (1978) are given in Appendix A3. 

For the present comparison, in addition to above properties for 
this track, following parameter has been assumed : 



Vertical Stress. (J v (kPa) 



Fig. 3.63 Comparison of results for vertical stress along 

sleeper at 18 cm depth for FAST Track. 
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Fig. 3.64 Comparison of results for rail displacement alon 

rail direction for FAST Track. 



100 MPa 


The range of parameters used for ELASTIC. F and 2D8N models is 
as given below : 


Tie Spacing = 0.508, 0.610, 0.711 m 
Ballast depth = 0.152, 0.305, 0.610 m 
Young's modulus of subgrade =13.79, 27.85, 55.16 MPa 

The parameters used for 3D20N program are given as below : 

Tie Spacing = 0.508, 0.711 m 
Ballast depth =0.305, 0.610 m 
Young's modulus of subgrade = 27.85, 55.16 MPa 

Figures 3.65 to 3.67 show the variation of maximum 
vertical subgrade pressure with tie spacing, ballast depth and 
Young's modulus of subgrade for all 4 models viz. So (1978) , 
ELASTIC. F, 2D8N and 3D20N models. 

Results predicted by 3D20N and ELASTIC. F model compare 
well with So (1978) model. 

3.5.3 VIENNA ARSENAL SINGLE SLEEPER AND DERBY PANEL FULL SCALE 
MODEL TESTS 

Janin et al. (1982), Janin et al. (1983a), Janin et al. 
(1983b) and Janin et al. (1983c) report the results of tests 
carried out at Vienna Arsenal and Derby. They also report the 
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Fig. 3.65 



Fig. 3.66 



Fig. 3.65-3.67 Comparison ot results withSO(1978) 






prediction of results by using elastic and elasto-plastic 
finite element model 'ROSALIE ' and a model based on layered 
elastic theory 'BISTRO '. 

Tests at Vienna Arsenal were carried out on five 
structures. For structures 1 to 4 , a soil type with E g = 20 MPa 
and for structure 5 , a soil type with E = 50 MPa were used. 
These tests were carried out in box type model and load was put 
on a single sleeper. 

Tests at Derby were carried out on four other different 
structures. All structures were tested for two different soil 
types viz. Boulder clay (E g = 11 MPa) and Keuper Marl (E s = 3 0 
MPa) . These tests were carried out on full scale panel with 
proper wheel loading. Structures 1 to 5 and structures 1 to 4 
tested at Vienna and Derby have been shown in Fig. 3.68 
respectively. Details of tracks used for Derby and Vienna tests 
are given in Appendix A4. 

Results from Vienna and Derby tests and their model 
predictions viz. elastic layered model 'BISTRO ' and finite 
element elastic and elasto-plastic 'ROSALIE ' models have been 
compared with predictions by models ELASTIC. F, 2D8N and 
3D20N. For ELASTIC. F model, for Derby track, sand layer has 
been assumed to have same Young's modulus as subgrade. Vienna 
single sleeper test could not be predicted by ELASTIC. F. Also 
Vienna tracks showing stiff layer over soft layer could not be 
predicted by this method. 

Figures 3.69 and 3.70 show theoretical and test results 
for vertical stresses and total resilient displacement for 
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Pig- 8.68 Structures tested at Vienna and Derby CJanin et al 1 983BD . 
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g. 3.69 Comparison of results for total resilient displace 
ment for Derby Panel tests 
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ig. 3.69 Comparison of results for total resilient displace 

ment for Derby Panel tests 
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Fig. 3.70 comparison of Results for Major principal 

stress for Deroy t-anel tests 


track panel tests at Derby with ELASTIC. F, 2D8N and 3D20N 
models. The results compare reasonably well. 

Figure 3.71 summarizes the test results in terms of 
vertical stress carried out by using single sleeper at Vienna 
and panel tests at Derby. Predictions, in general, fall within 
the data scatter except for Derby structure 1. 

Figure 3.72 compares horizontal stress distribution as 
measured for Vienna structure No. 2 and the same predicted 
inside the subgrade using 3D20N model. Figure 3.73 compares 
horizontal stresses on top of subgrade for Derby structures 
predicted by elasto-plastic 'ROSALIE ' analysis with 3D20N 
elastic analysis. 

Figure 3.74 shows the results obtained using 
elasto-plastic model 'ROSALIE ' for the purpose of preparation 
of design charts for Vienna type structures. However, due to 
insufficient data reported by Janin et al. (1983c) regarding 
proportion of depths of different layers, it has been decided 
to compare these results with the results of 3D20N for Vienna 
structures with full panel. 

From this comparison, it is clear that 3D20N model 
predicts vertical stresses throughout the track and horizontal 
stresses inside the subgrade well within the range with test 
results and with the analysis carried out by elasto-plastic 
analysis. Also, it is relevant to mention here that Desai and 
Siriwardane (1982) have recommended the use of two dimensional 
and three dimensional linear elastic finite element model for 
general use and have suggested that elasto-plastic 3-D finite 
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Fig. 3.71 Comparison of results in terms of vertical 

stresses for Vienna and Derby Tests . 
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Fig. 3.72 Comparison of results for horizontal stresses 

for Vienna structure NO. 2 
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Fig. 3.74 Comparison ot results with Rosalie elasto plastic 
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element analysis may be used in special circumstances such as 
evaluation of stresses on curves, loss of ties, higher axle 
loads etc. 

3.5.4. FIELD MEASUREMENTS BY RDSO ON INDIAN RAILWAYS 

Field measurements carried out by RDSO have been reported 
by Kapoor et al. (1967a), Kapoor et al. (1967b) and Mukerji and 
Yog (1974) . Measured results have been compared with the 
analytical results obtained by ELASTIC. F, 2D8N and 3D20N models 
in terms of deflection of sleeper, top formation pressure and 
the effect of double wheel load. In addition, field measurement 
of formation pressure due to the long haul train trials at 
Mandah road station, Allahabad division, 1986 have also been 
compared [RDSO (C) ] . 

Details of track used for tests carried out by RDSO along 
with the data assumed for purpose of modelling are given in 
Appendix A5. The results of these tests are also mentioned 
therein. 

3. 5. 4.1 Tests conducted on section No. 2 on Ambala - Ludhiana 
section (1964-1967) 

Comparison of results for section no. 2 on Ambala - 
Ludhiana section has been tabulated in Table 3.10. Pressure has 
been calculated at midheight in sand layer. 
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Table 3.10 Comparison of results for section no. 2, 
Ambala-Ludhiana section. 























3. 5. 4.2 Tests conducted at section No. 6 on Ambala Ludhiana 
section (1964-67) 

Comparison of results for section no. 6 on Ambala - 
Ludhiana section has been shown in Table 3.11. The predicted 
results seemingly do not compare well with the test results for 
this section. However, it is noted that range of scatter 
observed in this field study is very high and average value 
does not have much significance. 

3.5.4. 3 Tests conducted on a dead-end siding near canal bed 
site at RDSO, Lucknow 

Details of track structure for tests conducted at this 
site is also given in Table 3.12 in addition to Appendix A5. 
Comparison of results is shown in Table 3.13. 

3. 5. 4. 4 Field measurement of formation pressure due to long 
haul train trials at Mandah road station, Allahabad division 

Comparison of Results in terms of top formation pressure 
is given in Table 3.14. 

3.5.5 EMPIRICAL FORMULAE AVAILABLE IN LITERATURE 

Results in terms of top subgrade pressure from three model 
viz. ELASTIC. F, modified 2D8N and 3D2 0N models have been 
compared with the results obtained from empirical formulae 
available in literature [Kapoor et al. (1967a)] in Table 3.15 
for the track tested by RDSO with 0.20 m depth of formation for 



Table 3.11 Comparison of results for section no. 6, 
Ambala-Ludhiana section. 


Result in 

Observed 

Prediction by 

terms of 

Range 

Average 

elastic. f 

2D8N 

mod. 2D8N 

3D20N 

Formation 

Pressure 

(kPa) 

52.97 

to 

105.0 

84.37 

47.58 

25.01 

60.82 


Displacement 

(cm) 

0.736 

to 

0.152 

0.1016 



0.073 

— 

0.010 


Table 3.12 Track properties for canal bed site, Lucknow 
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Table 3.13 Comparison of results for canal bed site, Lucknow. 
(All pressure values are in kPa) 


Track 

structure 

type 

Average value 
of formation 
pressure observed 

Predicted by 

elastic. f 

mod. 2D8N 

3D20N 

1 

131.45 

109.87 

152.06 

114.80 

2 

63.76 

68.70 

58.86 

73.60 

3 

99.08 

95.16 

137.34 

101.00 

4 

55.92 

61.80 

56.90 

65.70 

5 

87.31 

82.40 

127.53 

87.31 

6 

49.05 

55.92 

55.92 

59.84 


Table 3.14 Comparison of results for Mandah Road Station, Allahabc 
division . (All pressure values are in kPa) . 


Speed 

(kmph) 

Average value 
of formation 
pressure observed 

Predicted by 

elastic. f 

mod. 2D8N 

3D20N 

65 

193.26 

137.34 

206.0 

116.7 


* Predicted for axle load = 22.8 t (for Box wagon) + 44.5 % dynamic 
augment for 66 kmph speed [RDSO (D) ] . 























Table 3.15 Comparison of results with empirical formulae 


Serial No. 

Empirical Method 

%t 

(kPa) 

1 

Observed 

86.32 

2 

Formula T.P. No. 245 

198.16 

3 

Formula Monograph no. 12 

53.96 


(average) 


4 

Formula Monograph no. 12 

313.92 


(maximum) 


5 

Talbot's Formula 

214.84 

6 

JNR formula (Japanese) 

149.11 

7 

German Formula 

206.01 

8 

Elastic. f 

69.65 

9 

Mod. 2D8N 

29.43 

10 

ID20N 

69.65 




section no. 2 on Ambala - Ludhiana section. The wheel load and 
track modulus used were 70 kN and 28.2 MPa respectively. 

From Table 3.15, it may be seen that results obtained by 
using empirical methods are 2 to 4 times greater than the 
measured results for this Indian track. While for the same 
track, ELASTIC. F and 3D20N model predict stresses within 20 % 
of the measured stresses which is the range of measurement for 
these results. 

3.6 COMPARISON OF THE THREE MODELS 

In this section, results obtained from 3 models viz. 
ELASTIC. F, 2D8N and 3D20N model for TRACK 1 have been 
discussed. 

It may be seen from Figs. 3.4 0 and 3.60 that results 

obtained from 2D8N lie within 10 % of results obtained by 3D20N 

model. But it may be noted that 2D8N will underpredict stresses 

on top of subgrade for smaller depth of formation and will 

overpredict it for large depths due to 2-D nature of this 

model. Also from Figs. 3.5 and 3.60, it may be noted that 

results obtained from ELASTIC. F lie from 5-20 % of results 

obtained by 3D20N m general and can vary up to 90 % for high 

values of E and low values of d, . 

s t 

The percentage variation in results predicted by each of 
the 3 models for same variation of different parameters has 
been compared in Table 3.16. Comparison of results for 
variation of different parameters viz., 


tie spacing, ballast 



Table 3 16 Comparison of results predicted by three 



COD 
































depth and E with results of So (1978) has already been shown 

5 

in Figs. 3.65 to 3.67. 

CPU time for the 3 models is compared in Table 3.17. 



Table 3.17 Comparison of CPU time on Convex 220 computer systems 
for three models 


Model type 

Elastic. f 

Mod. 2D8N 

3D20N 

CPU (secs) 

0.04 

1.6 

1774 





